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December  1986 

Chairman:  Dr.  Thom  J.  Hodgson 

Cochairman:  Dr.  Louis  Martin-Vega 

Major  Department:  Industrial  and  Systems  Engineering 

This  dissertation  addresses  the  problem  of  transporting  personnel, 
equipment  and  supplies  over  a long  distance  using  military  transport 
aircraft  (A/C),  of  which  one  or  more  may  need  to  use  aerial  refueling. 

It  is  desired  to  assist  the  decision  maker  in  answering  the  following 
questions:  (1)  How  many  A/C  are  needed?  (2)  How  should  the  total  load 

be  divided  among  all  the  used  A/C?  (3)  What  should  the  coordinates 
(location)  of  each  refueling  point  be?  (4)  What  should  the  amount  of 
initial  fuel  be  in  order  to  minimize  the  total  fuel  consumption  used  to 
complete  the  move? 

The  aim  of  this  dissertation  is  to  analyze  the  general  A/C  mid-air 
refueling  problem,  and  to  use  the  insights  obtained  from  the  analysis  to 
characterize  and  find  the  optimal  solution.  It  is  found  that  even 
though  a mathematical  model  for  the  general  problem  can  be  formulated, 
it  would  provide  little  insight  and  would  be  too  complicated  to  solve. 
Instead,  the  analysis  approach  used  in  this  dissertation  consists  of 
breaking  the  general  problem  down  into  smaller  and  simpler 
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(sub)problems  that  are  addressed  in  order  of  increasing  complexity. 
First,  the  various  relationships  between  fuel  consumption,  distance 
traveled,  cargo  weight,  and  initial  fuel  for  a typical  jet  A/ C are 
derived.  Next,  formulations  and  solution  procedures  of  the  first  five 
(sub) problems  are  presented.  These  include  both  the  problem  of  many  A/ C 
with  one  refueling  per  A/C  required,  and  the  problem  of  one  A/C  with  two 
refuelings  required.  The  spherical  convexity  of  these  problems  is 
shown. 
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CHAPTER  I 


INTRODUCTION 

Problem  Description  and  Model 

Consider  the  following  scenarios:  A decision  is  made  to  transport 

military  personnel,  equipment  and  supplies  from  an  air  force  base  in  the 
United  States  to  Cairo,  Egypt  to  aid  in  the  deployment  of  a squadron  of 
F-15  fighters  to  Egypt.  A similar  decision  could  have  been  made  to 
provide  needed  military  supplies  as  in  the  case  of  the  airlift  to  Israel 
during  the  1973  Arab-Israeli  War.  A third  case  might  be  a Rapid 
Deployment  Force  where  military  personnel  and  light  equipment  are 
transferred  to  a sensitive  area  to  intervene  in  case  of  an  emergency, 
e.g. , the  invasion  of  Grenada  in  1983. 

All  of  the  above  situations  are  similar  in  the  following  sense. 
Military  transport  aircraft  (A/C)  are  used  to  effect  the  move.  The 
equipment  and  supplies  are  either  packaged  on  pallets  or  consist  of 
rolling  equipment  large  enough  to  be  transported  as  self-contained 
units.  Equipment  and  men  are  loaded  into  the  transport  A/C,  taking  into 
consideration  various  physical  restrictions.  These  restrictions  include 
factors  such  as  type  and  availability  of  A/C,  as  well  as  loading 
considerations  such  as  maximum  weight  limit,  type  and  size  of  equipment, 
and  many  others. 

Given  a deployment  situation,  another  problem  arises  from  the  fact 
that  the  distance  between  the  origin  base  and  the  desired  destination  is 
sometimes  so  great  that  it  is  necessary  for  the  transport  A/ C to  refuel 
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en  route.  Many  times,  there  are  air  bases  along  the  way  where  the 
transport  A/C  could  land  and  refuel.  However,  refueling  may  be  per- 
formed more  efficiently  in  mid-air  using  special  A/C  called  "tankers," 
designed  for  this  purpose.  The  tanker  flies  to  a point  where  it  meets 
the  transport  A/C  and  supplies  it  with  fuel  while  both  are  in  the  air. 
The  objective  is  to  complete  the  move  so  as  to  minimize  the  total  fuel 
consumption  of  all  the  A/C  involved. 

In  analyzing  such  a problem,  the  following  decision  variables  need 
to  be  specified: 

1.  number  of  transport  A/C  of  each  type  that  are  needed 

2.  weight  of  the  load  going  into  each  A/C 

3.  initial  fuel  each  A/C  should  carry 

4.  transport  A/C  route 

5.  tanker  A/C  route 

Specifying  4 and  5 is  equivalent  to  determining  the  location  of  all 
refueling  points.  Figure  1-1  and  1-2  show  some  examples  of  possible 
routes . 

Figure  1-1  illustrates  the  case  of  two  tanker  bases.  The  transport 
flies  to  the  first  refueling  point  where  it  meets  a tanker  coming  from 
the  first  base  and  refuels  in  mid-air.  After  aerial  refueling  is  fin- 
ished, the  tanker  flies  to  its  base  while  the  transport  flies  to  the 
second  refueling  point  where  it  meets  another  tanker  from  the  second 
base.  Then,  after  refueling,  the  transport  heads  directly  for  the 
desired  destination.  Figure  1-2  depicts  a case  where  there  is  only  one 
tanker  base,  but  aerial  refueling  is  done  twice  using  tankers  origi- 


nating from  this  base 
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Tanker  Base  1 


0ri9in  Base  Destination 


0 Refueling  Point 

Transport  A/C  route 

— Tanker  route 


Figure  1-1: 


Illustration  of  aerial  refueling  and  possible  routes 
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Tanker  Base 


* 
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Destination 
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Transport  A/C  route 

Tanker  route 


Figure  1-2:  A second  example  of  a possible  route 
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The  general  problem  involves  many  variables,  such  as  different 
types  of  A/C,  the  possibility  of  many  tanker  bases,  and/or  multiple 
refueling,  etc.  Even  though  a mathematical  model  for  the  general 
problem  could  be  formulated,  it  would  provide  little  insight.  Rather, 
the  approach  taken  has  been  to  break  the  problem  down  into  smaller 
problems,  all  of  which  are  special  cases  of  the  general  problem.  By 
starting  with  the  simplest  problem  and  modelling  them  in  increasing 
order  of  difficulty,  one  captures  the  essence  of  the  general  problem. 

The  simpler  problems  also  serve  as  building  blocks  for  subsequent,  more 
complicated  ones.  Figure  1-3  illustrates  the  various  problems  which  are 
described  below. 

Problem  (Pi).  There  are  N A/C,  and  any  A/ C used  is  required  to  fly 
directly  to  the  destination.  It  is  desired  to  determine  how  many  A/ C to 
use  and  how  to  divide  the  load  among  them. 

Problem  (P2).  There  is  only  one  A/C  with  predetermined  cargo 
weight,  but  it  is  unable  to  fly  all  the  way  to  the  destination  without 
refueling.  The  location  of  the  optimal  refueling  point  and  the  initial 
fuel  required  for  both  the  transport  A/C  and  the  tanker  A/ C need  to  be 
determined. 

The  first  two  problems  are  fundamental  because  no  other  problem  can 
be  solved  without  knowing  first  how  to  solve  these  two. 

Problem  (P3).  This  problem  has  the  complexity  of  (PI)  and  (P2) 
combined.  Here  there  are  two  A/C.  The  cargo  weight  that  goes  into  each 
one  needs  to  be  determined  and,  depending  on  that  weight,  the  A/C  might 
need  to  refuel.  Therefore,  the  location  of  the  optimal  refueling  point 
of  each  A/C  (if  refueling  is  necessary)  needs  to  be  specified  as  well. 
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(PI) 

N A/C 

No  refueling 
required 

(1  leg) 


Figure  1-3: 


Problems  and  Solution  Stages 


(P2) 


(P5) 


Different  cases  of  the  general  problem  and  the  solution 
stages. 


7 


Problem  (P4).  This  problem  is  essentially  the  same  as  problem  (P3) 
with  N A/C  instead  of  two. 

Problem  (P5).  This  problem  is  an  extension  of  (P2)  to  include  the 
possibility  of  refueling  twice  instead  of  once. 

Problem  (P6).  This  problem  can  be  considered  as  an  extension  of 
(P5)  to  include  two  A/ C instead  of  one  or  as  an  extension  of  (P3)  with 
two  refuelings  instead  of  one. 

Problem  (P7).  This  problem  is  similar  to  (P6)  but  with  N A/C.  If 
this  problem  is  solved,  then  the  extension  to  the  case  of  multiple 
refuelings  should  be  easy  and  obvious. 

In  Figure  1-3,  the  arrows  indicate  the  sequence  that  must  be 
followed.  In  order  to  solve  a certain  problem,  one  must  first  solve  all 
the  problems  leading  to  it. 


Literature  Survey 

There  exists  no  research  in  the  reported  literature  that  deals  with 
this  problem  in  its  generality,  i.e.,  (P7).  However,  by  breaking  it 
down  into  subproblems,  one  can  identify  a fundamental  and  important 
subproblem  which  involves  one  transport  A/C  where  aerial  refueling  is 
required  (i.e.,  P2).  This  subproblem  concerns  itself  with  finding  the 
location  of  the  refueling  point  such  that  the  total  fuel  consumption  is 
minimized  [i.e.,  a location  problem  on  a sphere  (earth)].  The  objective 
function  is  nonlinear  in  the  distance,  and  the  constraints  (the  A/C 
cannot  be  allowed  to  run  out  of  fuel)  depend  on  some  of  the  decision 
variables. 

A closely  related  problem  is  found  in  the  literature  of  location 
theory,  but  it  has  a linear  objective  function.  This  problem  is  the 
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spherical  Weber  problem  with  maximum  distance  constraints.  In  the  Weber 
problem,  it  is  necessary  to  determine  the  optimal  location  of  a source 
so  as  to  minimize  the  sum  of  weighted  distance  (minisum  objective)  from 
the  source  to  a finite  number  of  given  demand  points  whose  location  is 
known. 

The  spherical  Weber  problem  on  a unit  sphere  (without  maximum 
distance  constraints)  can  be  stated  mathematically  as  follows: 

m 

Minimize  f(X)  = 2 B • d(X,P.)  (1-1) 

i=l  1 1 

S.T.  | | X | |2  = 1,  and  X e E3  (1-2) 

where 

f*  Vi 

P^  = location  of  i1"  demand  point  on  the  surface  of  the  unit 

sphere,  i = 1 , . . . , m 

3^  = weight  corresponding  to  the  i1"*1  demand  point, 

3^  > 0 , i = 1 , . . .,m 

X = location  of  the  source 

d(X,P^)  = distance  between  the  source  and  the  itlri  demand  point, 

measured  along  the  surface  of  the  sphere 

| |*  | | = the  Euclidean  norm. 

Several  researchers  have  investigated  the  nature  of  the  optimal 
solution  for  the  unconstrained  spherical  Weber  problem  (a  summary  of 
this  can  be  found  in  Appendix  C).  Aly  et  al.  [1]  showed  that  the  search 
for  an  optimal  solution  to  this  problem,  where  demand  points  are  not 
located  entirely  on  a great  circle  arc  (see  Appendix  B for  the 
definition  of  a great  circle),  can  be  restricted  to  the  spherically 
convex  hull  of  the  demand  points.  Drezner  [10]  proved  that  if  the 
demand  points  are  located  on  a great  circle,  so  is  the  optimal  solution 
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point.  Drezner  and  Wesolowsky  [12]  gave  properties  and  results  for  the 
problem,  some  of  which  were  shown  to  be  invalid  by  Ayken  [3].  He 
restated  those  properties  and  results  and  established  additional  ones. 
One  of  the  restated  results  is  that  if  all  the  demand  points  of  a 
spherical  Weber  problem  are  included  within  a spherical  disk  of  radius 
less  than  or  equal  to  n/4,  every  minimum  of  the  problem  is  a global 
minimum.  An  additional  result  was  as  follows:  Suppose  g^,  i = 

1,  . . .,  m are  the  given  points  and  their  respective  weights,  where  m _>. 
3 and  3^  _>.  0*  If  all  the  demand  points  are  included  within  a spherical 
disk  of  radius  tt/2,  and  at  least  three  demand  points  are  not  spherically 
colinear,  the  spherical  Weber  problem  possesses  a unique  minimum. 
Furthermore,  that  minimum  is  in  the  spherical  disk  covering  the  demand 
points. 

Katz  and  Cooper  [20]  gave  a condition  under  which  a demand  point 
(9  ,<J>  ) is  a local  minimum.  Let  a.,  represent  the  shortest  arc  distance 

1C  tC  IK. 

between  P^  and  P^,  where 

cos  aik  = cos  ^ cos  <j>k  cos(ek  - Q ) + sin  ^ sin  (1-3) 

The  demand  point  P^  is  a local  minimum  if 


(1-4) 


m g 


i*k 


+ cos  <j>,  sin  a] 
k i 


(1-5) 


Z 

i=l 

i*k 


m ^ 


m 


(1-6) 


3,  = the  weight  associated  with  P . 

K.  K. 
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This  condition  is  similar  to  the  so-called  Kuhn  modified  gradient 
condition  for  the  planar  case  (see  Francis  and  White  [16]). 

Several  solution  procedures  to  the  spherical  Weber  problem  without 
maximum  distance  constraints  have  been  developed.  Litwhiler  and  Aly 
[22]  solved  the  problem  using  two  methods.  The  first,  the  Map  Projec- 
tion Algorithm  (MPA),  projects  the  current  solution  point  to  the  plane 
using  the  azimuthal  equidistance  projection,  solves  the  problem  in  the 
plane  using  euclidean  distance,  then  transforms  it  back  to  the  sphere. 
This  is  slow  due  to  the  back  and  forth  transformation.  The  other 
method,  Cyclic  Search  Algorithm  (CSA),  searches  in  a cycle  in  a direc- 
tion orthogonal  to  a meridian  along  a great  circle  track  and  then  along 
a meridian.  However,  the  North  Pole  has  to  be  transformed  to  the  cur- 
rent solution  point  in  order  to  find  the  great  circle  track.  This 
algorithm  is  derivative  free  and  is  faster  than  the  MPA. 

Katz  and  Cooper  [20]  and  Drezner  and  Wesolowsky  [12]  developed 
similar  procedures  for  the  Weber  problem.  Both  are  iterative  procedures 
that  use  the  derivative  of  the  objective  function.  The  difference 
between  the  two  algorithms  is  in  the  functional  form  used  to  char- 
acterize the  great  circle  distance. 

Ayken  [3]  also  developed  two  algorithms.  The  first  is  a derivative- 
free  algorithm  called  Cyclic  Meridian  Parallel  Search  (CMPS).  It  is 
similar  to  the  CSA  of  Litwhiler,  but  faster  because  no  transformation  of 
poles  to  the  current  point  is  needed.  Also,  CMPS  can  be  used  to 
minimize  any  type  of  function  restricted  to  any  surface  on  which  the 
function  and  the  point  can  be  easily  represented  by  curvilinear  coordi- 
nates. The  second  algorithm  is  called  the  Geodesic  Descent  Algorithm 
(GDA).  The  search  direction  of  GDA  is  the  geodesic  curve  whose  tangent 
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vector  is  the  gradient  direction  of  the  objective  function  projected 
onto  the  tangent  plane.  GDA  is  analogous  to  the  steepest  descent  algo- 
rithm used  for  the  unconstrained  problem  in  En,  and  GDA  is  faster  than 
either  CMS  or  MPA.  Also,  like  CMS,  GDA  can  be  used  to  optimize  any 
function  defined  on  any  surface  for  which  the  following  conditions  are 
met:  it  is  possible  to  project  a point  lying  outside  the  surface  onto 

it;  it  is  straightforward  to  find  the  path  of  the  geodesic  descent;  and 
partial  derivatives  of  the  function  are  easily  found.  Note  that  the 
iterative  scheme  proposed  by  Katz  and  Cooper  [20]  could  be  considered  as 
a fixed-step  version  of  GDA  as  shown  by  Ayken  [3]. 

All  of  the  above  algorithms  solve  the  spherical  Weber  problem  with- 
out maximum  distance  constraints.  They  converge  to  a local  optima 
(which  is  also  a global  optima)  if  the  search  region  is  limited  to  a 
spherical  disk  of  tt/4  radius.  Moreover,  the  optimal  point  may  not  be 
unique. 

Now,  we  turn  our  attention  to  the  constrained  version  of  the  pro- 
blem. Constrained  location  problems  in  En  have  been  studied  by  several 
researchers.  Love  [23]  proposed  a scheme  to  handle  various  types  of 
spatial  constraints  which  make  use  of  convex  programming  and  penalty 
function  techniques.  Love  and  Morris  [24]  studied  the  computational 
aspects  of  the  solution  of  constrained  multifacility  location  problems 
involving  Lp  distances  using  nonlinear  optimization  techniques.  Hurter 
et  al.  [18]  provided  some  properties  of  constrained  location  problems 
when  the  distance  is  derived  from  a norm. 

Another  version  of  the  constrained  location  problem  is  one  that 
includes  the  presence  of  a forbidden  region  in  which  no  path  is  per- 
mitted to  enter.  This  was  studied  by  Katz  and  Cooper  [6,19,21]  using  Lp 


distance 
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Ayken  [3]  formulated  the  spherical  Weber  problem  with  maximum 
distance  constraints  in  which  the  location  of  the  sources  is  kept  within 
a distance  St  of  a fixed  point  Ut,  where  Ut  may  be  one  of  the  demand 
points.  The  maximum  distance  constraint  on  a unit  sphere  can  be  char- 
acterized as  follows: 

Let  I = {X  : | |X | | 2 = 1 , Arc  cos  (X'l^)  < S , Xe  E3} 

Inequality  arc  cos  (X'Ut)  <_  St  is  equivalent  to  the  following  two 


inequalities : 

-X'ut  <_  - cos  (St) 

(1-7) 

X'ut  1 1 

(1-8) 

The  second  inequality  is  redundant  so  the  spherical  Weber  problem  on  a 
unit  sphere  with  maximum  distance  constraints  is  formulated  as  follows: 


m 

Min  f(X)  = E 0 • d(X,P.)  (1-9) 

i=l  1 1 

S.T.  | | X | |2  = 1 (1-10) 

-X'ut  <_  -Dt  t = 1,  . . . , N (1-11) 

XeE3  (1-12) 

where  Dt  = cos  (St). 


Note  that  this  formulation  has  a linear  objective  function  and  all 
points  are  expressed  using  cartesian  coordinates  rather  than  spherical 
coordinates.  Ayken  [3]  solved  this  problem  by  first  finding  a vector 
tangent  to  the  improving  feasible  direction.  This  is  done  by  solving 
the  following  linear  program  subproblem. 
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Let  J = It  : U'  X = D I 
1 t tJ 


Minimize  f(x)' 

• d 

(1-13) 

S.T.  -U'  d < -D 
t - t 

t e J 

(1-14) 

X'd  = 0 

(1-15) 

-1  < d.  < 1 

i = 1,2,3 

(1-16) 

where  di  = the  ith 

component  of  d. 

The  component  d^  can  be  solved  in 

terms  of  the  other  components  of  the  vector  d (i.e.,  and  d2)  using 
equation  (1-15).  This  reduces  the  subproblem  and  makes  it  possible  to 
solve  graphically.  Before  conducting  a line  search  along  the  direction 
d,  the  maximum  step  size  that  can  be  taken  from  the  current  solution  in 
this  direction,  without  violating  any  constraint,  has  to  be  found.  A 
formula  is  given  for  getting  the  maximum  step  size.  Then  a line  search 
in  this  direction  is  conducted  to  find  a new  point.  The  process  is 
repeated  until  convergence  is  reached. 

Another  version  is  the  constrained  maxmin  or  minmax  facility 
location  problem  on  the  sphere.  Several  researchers  have  worked  on  this 
problem,  including  Drezner  and  Wesolowsky  [13,14]  and  Drezner  [11]. 

Darnell  and  Loflin  [7,8]  and  Waite  [28]  attacked  a simplified 
version  of  problem  (P2)  and  (P5)  in  which  there  are  unrestricted  but 
finite  number  of  fixed  refueling  points  with  different  fuel  costs.  The 
objective  is  to  fly  the  predetermined  tour  with  minimum  fuel  cost. 
Darnell  and  Loflin  [7,8]  used  linear  programming  to  solve  the 
deterministic  version  of  this  problem.  While  Waite  [28]  concentrated  on 
the  stochastic  nature  of  this  problem  and  solved  it  using  stochastic 
dynamic  programming. 
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Assumptions 

In  order  to  focus  on  the  analytical  elements  of  this  study,  the 

following  simplifying  assumptions  are  made. 

1.  The  total  cargo  to  be  transported  can  be  divided  in  any  way  desired, 
i.e.,  one  can  ignore  the  fact  that  the  cargo  is  packaged  in  pallets 
and  that  one  does  not  have  the  freedom  to  split  the  cargo  load  in 
any  way  except  as  the  weight  of  the  individual  pallets  permits.  In 
essence,  it  is  assumed  that  the  total  cargo  weight  is  a continuous 
variable.  The  weight  associated  with  the  assignment  of  people  to  an 
aircraft  will  also  be  assumed  to  be  a continuous  variable. 

2.  Only  the  following  limitation  on  the  weight  is  to  be  considered: 
the  gross  weight  of  the  A/C  should  not  exceed  the  maximum  take-off 
weight  of  the  A/C  (MTOW),  [i.e.,  A/C  empty  weight  + cargo  weight  + 
fuel  weight  <_  MTOW]. 

3.  There  are  no  load  balance  or  size  restrictions.  That  is,  any  part 
of  the  cargo  can  fit  into  any  part  of  the  A/C  as  long  as  the  maximum 
take-off  weight  restriction  is  not  violated.  This  is  reasonable  in 
most  cases  because  almost  all  pallets  can  fit  anywhere  in  the 
transport  A/C  except  possibly  in  the  rear-most  position  of  the  A/C. 

4.  The  earth  is  a perfect  sphere  with  radius  R,  where 
R = 3,920  statute  miles 

= 3,404  nautical  miles. 

5.  The  A/C  follows  the  great  circle  route  in  flying  from  one  point  to 
another. 

6.  Weather  at  altitude  (i.e.,  the  jet  stream)  is  assumed  to  be 
negligible.  Inclusion  of  this  factor  would  be  relatively 
straightforward  but  would  unnecessarily  complicate  the  analytic 


development 
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7.  Aerial  refueling  takes  a negligible  amount  of  time.  Thus,  the 
region  in  which  fuel  transfer  takes  place  is  considered  to  be  a 
single  point. 

8.  Everything  takes  place  at  altitude.  Thus,  fuel  consumed  for  take- 
off is  ignored. 

9.  Only  one  type  of  transport  A/C  is  available. 

Scope  of  the  Study 

This  study  is  concerned  with  analyzing  the  general  A/ C mid-air 
refueling  problem  described  earlier  in  order  to  characterize  and  find 
its  optimal  solution.  The  limitations  of  this  study  are  discussed  in 
the  previous  section  under  "Assumptions." 

The  analysis  approach  used  in  this  dissertation  consists  of 
breaking  the  general  problem  down  into  smaller  and  simpler  problems  that 
are  addressed  in  order  of  increasing  complexity  up  to  the  most  general 
one.  Starting  with  problems  (PI)  and  (P2),  the  first  five  problems  are 
solved  in  the  sequence  indicated  by  the  arrows  in  Figure  1-3. 

Chapter  II  contains  derivation  and  characterization  of  the  fuel 
consumption  function  and  other  related  functions.  These  functions 
describe  the  relationship  between  fuel  consumption,  distance  traveled, 
cargo  weight,  and  initial  fuel  of  a typical  jet  aircraft. 

Chapter  III  contains  a formulation  of  problem  (PI).  An  exact 
solution  for  this  problem,  showing  the  optimal  number  of  A/ C to  use, 
their  optimal  cargo  weights  and  their  initial  fuel,  is  provided. 

Chapter  IV  contains  an  extensive  description  and  formulation  of 
problem  (P2).  This  problem  requires  the  location  of  the  optimal 
refueling  point  for  one  A/C  to  be  determined  together  with  the  initial 
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fuel  for  both  the  transport  and  the  tanker  A/C.  Spherical  convexity  of 
this  problem  is  discussed  and  proved.  A procedure  to  find  the  optimal 
solution  is  also  provided. 

Chapters  V and  VI  present  formulation  and  solution  procedures  for 
problems  (P3)  and  (P4),  respectively. 

Chapter  VII  addresses  (P5).  A solution  procedure  based  on  the 
developments  of  Chapter  IV  is  proposed  and  discussed.  It  is  shown  to 
converge  at  least  locally  and  that  it  is  stable.  Computational  results 
of  some  test  problems  support  the  conjecture  that  this  solution 
procedure  converges  globally  for  real  world  problems.  Finally,  Chapter 
VIII  summarizes  the  results  obtained  in  this  dissertation  and  provide 


directions  for  future  research 


CHAPTER  II 


THE  FUEL  CONSUMPTION  FUNCTION  AND  OTHER  RELATED  FUNCTIONS 

Derivation  of  the  Fuel  Consumption  Function  and 
Other  Related  Functions 

Here,  mathematical  representation  of  the  various  relationships 
between  fuel  consumption,  distance  traveled,  cargo  weight,  and  initial 
fuel  for  a typical  jet  A/C  are  derived. 

The  United  States  Air  Force  [29]  provided  raw  data  describing  jet 
A/C  performance  at  different  gross  weights,  different  altitudes,  and 
different  speeds.  These  raw  data  have  to  be  transformed  before  they  are 
useful.  A sample  of  these  data  is  shown  in  Figure  2-1  for  a C-5A 
transport  A/C  flying  at  an  altitude  of  31,000  feet.  Different  curves 
are  presented  for  different  A/C  weights,  and  they  show  the  distance 
traveled  per  1 ,000  pounds  of  fuel  burned  at  a given  speed  and  at  that 
given  gross  weight  (GW).  Consider  the  points  on  the  curve  indicating  99 
percent  maximum  specific  range.  It  is  assumed,  without  loss  of 
generality,  that  the  transport  A/C  will  operate  at  those  points  since 
speed  is  maximized,  while  fuel  consumption  is  only  1 percent  greater 
than  the  minimum  possible. 

The  distance  traveled  in  miles  per  1,000  pounds  of  fuel  burned  when 
the  A/C  gross  weight  is  GW,  denoted  here  as  MPF  (GW),  is  plotted  against 
gross  weight  GW  in  Figure  2-2  using  the  99  percent  line  in  Figure  2-1. 
Least  square  linear  and  quadratic  fits  to  the  data  were  tried,  and  both 
fit  well  as  seen  in  Figure  2-2. 
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Air  Nautical  Miles  per  1,000  Pounds  of  Fuel 
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Model  C-5A 


SPECIFIC  RANGE 
4 Engines  31,000  Feet 


Figure  2-1:  Data  for  C-5A  aircraft 
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MPF(GW)  - Miles  per  1000  lbs.  of  fuel  burned  for 
a given  gross  weight 


Mi 1 1 i 1 1 1 ► 

320  400  480  560  640  720 

GW  = Gross  weight  in  1000  lbs. 


Figure  2-2:  The  MPF  function  versus  GW. 
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For  the  linear  fit, 

MPF  (GW)  = aQ  + a1  GW 
where 

a = 36.2829 
o 

ax  = -0.027 

and  p = correlation  coefficient  = -0.9919. 

For  the  quadratic  fit, 

MPF  (GW)  = bQ  + bL  GW  + b2GW2 

where 

b = 43.7616 
o 

bl  = -0.0576 

b2  = 2.94  x 10“5 

and  p = -0.9983. 

The  MPF  function,  in  both  forms,  will  be  used  to  find  the  distance 
the  A/C  can  travel  when  it  has  initial  fuel  gQ  and  cargo  weight  wQ. 

Also,  it  will  be  used  to  find  out  how  much  fuel  is  consumed  when  the  A/ C 
travels  a given  distance  d,  its  cargo  weight  is  wQ,  and  the  initial  fuel 

is  gQ. 

Note  that  at  any  moment  during  the  flight, 

GW  = EW  + wQ  + f (2-3) 

where  f = the  present  (instantaneous)  amount  of  fuel  and  EW  = A/ C empty 
weight.  During  flight,  f changes  due  to  fuel  burning;  thus,  GW  changes 
while  EW  and  wq  stay  the  same.  So, 


(2-1) 


(2-2) 


d GW  = d f 


(2-4) 
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The  Range  Function 

Let  R(go,wQ)  = the  range  of  the  tanker  A/C  when  Its  initial 
fuel  is  gQ  and  cargo  weight  is  wq. 

Initial  Gross  Weight  = EW  + wq  + gQ 
Final  Gross  Weight  = EW  + wQ 


then. 


EW  + w + g 

R(gQ,w  ) = f MPF  (GW)  d GW 

EW  + w 

o 


= r MPF  (GW)  d f 
0 


= f MPF  (EW  + w + f)  d f 
0 

(a)  Using  the  linear  fit  for  MPF,  this  results  in 
R(g0,wo)  = (a;  + alWQ  +4"  g0)  go 


where  a'  = a + a EW 
o o 1 


(b)  For  a quadratic  fit  of  MPF,  this  results  in 


(2-5) 


(2-6) 

(2-7) 


R^go,w0)  = [b; + biw0  + Vo 

bl  go 

+ (b2wo  + — >go  + b2  ~T]  go 

where  b*  = b + b,  EW  + b0EW2 
o o 1 2 

b;  = bl  + 2b2EW 


(2-8) 

(2-9) 

(2-10) 


Figure  2-3  shows  R(g0»wQ)  for  values  of  wQ  and  all  values  of  gQ 
using  both  linear  and  quadratic  fits. 
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g0=  Initial  fuel  in  1000  lbs. 


Figure  2-3:  Range  as  a function  of  initial  fuel  and  cargo  weight. 
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The  Fuel  Consumption  Function 

Now,  it  is  necessary  to  find  the  fuel  consumed  when  the  A/C  flies  a 
distance  d [d  <_  R(g0,wQ)]  with  initial  fuel  gQ  and  cargo  weight  wQ. 


Let 


Sf 

FC(go’wo’d) 


= final  amount  of  fuel  left,  and 

= fuel  consumed  when  initial  fuel  is  gQ , cargo  weight 
is  wQ,  and  the  distance  flown  is  d, 


then, 


FC(g0,w0,d)  = g0  - gf 


(2-11) 


Thus,  finding  FC  reduces  to  finding  gf,  which  can  be  accomplished  by 
solving  for  gf  in  the  following  equation: 


So 

d = / MPF  (EW  + w + f)  df 

o 

§f 

(a)  For  a linear  fit  of  MPF,  one  has 


(2-12) 


d = / [a  + a. (EW  + f + w )]  df 

o 1 o 

gf 

= ao(gQ  - gf)  + 3i(EW  + wo)(gQ  " gf)  + 4"  " gf} 


Solving  for  gf  results  in 


§f  = " 31  * 


V(a  + aiSQ)2  " 2aid 


a. 


(2-13) 
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Only  the  plus  sign  makes  gf  physically  possible,  so 


(2-14) 


Since  the  fuel  consumed  = FC(gQ,wo,d)  = gQ  - gf,  then 


(2-15) 


where 


a = a 


o 


+ ax(EW  + wQ) 


(2-16) 


Figure  2-4  shows  a typical  FC(go,wQ,d)  function  for  the  linear 

case. 

(b)  The  same  procedure  can  be  used  to  find  FC(g0,w0,d)  for  the  case  of 
a quadratic  fit  of  MPF.  The  solution  involves  finding  the  root  of  a 
cubic  equation  by  solving  for  gQ  in  the  following  equation: 


2 3 2 2 

3"  Sf+  — Sf  + ClSf 


(2-17) 


and  setting  FC(g0,w0,d)  = gQ  - gf 


where 


C1  = b0  + (EW  + w0)  + b2  (EW  + wQ)2 


(2-18) 


and 


C2  = 2 b2(EW  + wQ)  + b1 


(2-19) 
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d = distance  in  nautical  miles 


Figure  2-4:  Fuel  consumption  as  a function  of  distance  for  the  linear 

case 
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Using  the  CRC  Tables  [27],  the  solution  for  of  equation  (2-17) 
is  found  by  letting 


P 

q 

r 

el 

e2 

e3 

A 

B 


1 f2 

2 b2 

b2 


"3  Jo  , 2 2 . 

b7  {cigo  ~T  go 


1 2n 

T q ■ p ) 


-|y(2p3  - 9pq  + 27r) 


V 


e?  jh 

4 27 


\T 


2 + e3 


\r 


2 " ®3 


Then, 

gf  = A + B - £ 

The  Fuel  Requirement  Function 

Let  FN(wQ,d)  = the  exact  amount  of  fuel  needed  to  fly  a distance  d 
when  the  cargo  weight  is  wQ.  It  is  needed  to  find  FN. 

(a)  For  a linear  fit  of  MPF,  if  the  distance  d to  be  flown  and  the 
cargo  weight  wQ  are  known  in  advance,  then  the  amount  of  initial  fuel  gQ 
must  be  at  least  as  big  as  the  amount  of  fuel  needed.  If  gQ  is  set  to 
be  equal  to  FN(wQ,d),  then  we  can  use  the  range  function  R(g0,wQ)  to 
solve  for  gQ. 


That  is,  set 
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R(g0,w0)  = d = (a;  + aiwQ  + — gQ)  gQ 


(2-6) 


Solving  for  gQ  results  in  the  following  equation: 


FN(w0,d)  = gQ  = -wQ  - — + 


\j (a'  + a.w  + 2a.d 
y o 1 o 1 


(2-20) 


‘1  “1 

(b)  For  the  quadratic  fit  of  MPF,  the  following  equation  is  solved  for 
gQ  and  then  FN(wQ,d)  is  set  = gQ: 


Vo  + (Vo  + 4_)  So 


+ ( b ' + b.'  w + b„w  ) g - d = 0. 
o 1 o 2 o o 


(2-21) 


Again,  using  the  CRC  Tables  [27],  the  solution  for  gQ  of  equation 
(2-21)  is  found  by  letting 

bl 

P =W°+^ 

q = r-  (b'  + b! w + b„w3) 
b^  o 1 o 2 o 

. -d 


el  = T ^3q  “ p2) 


e2  = -|y(2p3  - 9pq  + 27r) 


/2  3 

® 1 

4 27 


2 ‘ e3 


2 e3 
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Then, 

gQ  = A + B - 

Figures  2-5  and  2-6  show  plots  of  FN  (wQ,d)  for  the  case  of  linear 
fit  of  MPF  as  a function  of  d and  wQ,  respectively.  Table  2-1 
summarizes  the  various  fuel  functions. 


Characteristics  of  the  Fuel  Functions 
The  Fuel  Requirement  Function;  FN  (w,d) 


FN  (w,d)  = -w  + 

h 


V< 


a'  + aw)  + 2a, d 
o 1 1 


(1)  The  function  FN  (w,d)  is  a strictly  convex  and  increasing  function 
of  the  distance  d for  any  given  value  of  the  cargo  weight  w. 

Proof : 

Take  the  first  and  second  derivatives  of  FN  with  respect  to  d and 
note  that  w here  is  fixed  at  a constant  value. 


3 FN 

3“T 


1 


+ a^)2 


+ 2a^d 


> 0 


3 d2  [ (a^  + a^)2  + 2a1d]3/2 

because  a^  < 0.  So,  FN  is  strictly  convex  in  d.  Moreover,  since  the 
first  derivative  is  always  positive,  then  FN  is  an  increasing  function 
of  d. 

(2)  The  function  FN  is  a strictly  convex  and  increasing  function  of  the 
cargo  weight  w for  any  given  value  of  d.  This  can  be  shown  in  a 
straightforward  fashion  as  in  the  previous  case. 
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d = distance  in  nautical  miles 


Figure  2-5:  Fuel  requirements  as  a function  of  distance  and 

cargo  weight. 
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w = cargo  weight  in  1000  lbs. 


Figure  2-6:  Fuel  requirements  as  a function  of  cargo  weight 

and  distance. 


TABLE  2-1 : Summary  of  the  fuel-related  functions 
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We  prove,  in  Appendix  D (Theorem  D-l),  that  if  the  distance  d is 
measured  along  the  surface  of  a sphere  (earth),  then  the  function  FN  (w,d) 
is  spherically  convex  (s-convex)  over  a spherical  disc  of  radius  < irM. 

The  Fuel  Consumption  Function:  FC  (g,w,d) 


(1)  The  function  FC  (g,w,d)  is  a concave  and  increasing  function  of  the 
distance  d for  any  given  values  of  the  initial  fuel  g and  the  cargo 
weight  w. 

(2)  The  function  FC  (g,w,d)  is  a convex  and  increasing  function  of  w 
for  any  given  values  of  g and  d. 

Again,  both  cases  can  be  shown  to  be  true  in  a fashion  similar  to 


FC(g,w,d)  = g + w + — - 

al 


a 


1 


that  of  FN(w,d) 


CHAPTER  III 


PROBLEMS  (PI)— MODEL  AND  SOLUTION 


Description  and  Formulation  of  (PI) 

A total  cargo  weight  of  W is  to  be  transported  from  one  base  to 
another  where  the  distance  between  them  is  Dqq»  There  are  N identical 
transport  A/C  available  at  the  origin  base  "0."  It  is  required  that 
each  of  the  M _<  N A/C  used  fly  directly  to  the  destination  "D." 

The  following  variables  need  to  be  determined: 

1.  How  many  transport  A/C  to  use,  i.e.,  M = ? 

2.  How  much  cargo  should  be  loaded  into  each  A/C,  i.e.,  wR  = ? for  n 
1 , 2 , . . . , M. 

3.  How  much  fuel  should  be  put  in  each  A/C,  i.e.,  gR  = ? for  n = 1,  2 

3 , . . . , M. 

The  objective  is  to  minimize  the  total  fuel  consumption. 

The  following  definitions  are  necessary: 

EW  = Transport  A/ C Empty  Weight 

MTOW  = Maximum  Take-off  Weight  of  the  transport  A/C 

Fmax  = Maximum  fuel  capacity  of  the  transport  A/C 

Rn  (gw)  = Range  of  the  nc^  transport  A/C  when  its  initial 

fuel  is  gn  and  cargo  weight  is  wn 
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FN  (w  ,d) 
n n’ 


= Exact  amount  of  fuel  needed  by  the  n*"^1  transport 
A/ C to  fly  a distance  d when  its  cargo  weight  is 


w_ 


Note  that  if  wQ  = 0 for  some  n,  then  that  A/C  will  not  be  used. 
Hence,  M can  be  determined  easily  once  wn's  are  known.  Also,  if  wn  > 0 
is  known,  then  the  amount  of  fuel  needed  (FNCw^Dq^))  can  be  determined 
easily.  Since  all  M <_  N A/C  used  are  flying  directly,  we  could  set  the 
amount  of  initial  fuel  gn  to  be  equal  to  the  amount  of  fuel  needed, 
i.e. , 


g 


n 


^wn,EW 


if  w > 0 
n 

if  wn  = 0 


N 

Thus,  total  fuel  consumption  would  be  equal  to  E g . Stating  the 

i n 

n=l 

problem  mathematically, 


N 


(PI)  Minimize  E 

n=l 

gn 

(3-1) 

N 

S.T.  E w = W 

i n 

n=l 

(3-2) 

gn  + wn  _<  (MTOW 

- EW) 

V n 

(3-3) 

0 < g < F 
— Bn  — max 

V n 

(3-4) 

Sn  = 


(wn,EW 


if  wn  > 0 


V n (3-5) 


if  wn  = 0 


wn>° 


V n 


(3-6) 
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This  problem  is  nonlinear  because  equation  (3-5)  is  nonlinear. 


Solution  of  (PI) 


The  following  two  theorems  tell  us  the  optimal  number  of  A/C  to  use 
and  their  optimal  cargo  weights. 


Theorem  3-1 

Given  that  M A/C  are  to  be  used,  the  total  cargo  load  should  be 
divided  equally  among  those  M A/C,  i.e.. 


n = 1,  2,  . . .,  M 
n = M + 1 , . . . , N 


Proof  of  Theorem  3-1.  Since  M A/C  will  be  used,  then 


wn  p 0 for  n = 1 , 2 , 
and  wn  = 0 for  n = M + 1 
Then  from  equation  (3-5),  it  must  be  true  that 


, M 
.,  N 


gn  FNn^wn,DOD^  for  n - 1>  2>  • • •,  M 

Thus,  the  problem  in  this  case  is  stated  as 


M M 


(Tl) 

Min  E gn  = 

n=l 

E 

n-1 

FN  (w  ,D) 
n n’  OD 

M 

S.T. 

y w = W 
n-1  n 

wn  > 0 

n = 

1,  • • 

.,  M 

Sn  + wn  £ (“TOW 

- EW) 

n = 

1 , • • 

.,  M 

(3-3) 

0 < g < F 
sn  -1  rmax 

n = 

1,  • • 

M 

(3-4) 

The  last  two  equations  will  be  ignored  for  the  moment.  If  the 
optimal  solution  to  the  resulting  relaxed  problem  satisfies  them,  the 
solution  is  optimal  to  the  original  problem  (Tl).  The  relaxed  problem 
(RT1)  is  as  follows: 

M 

(RT1)  Min  E FNn  (w^D^)  (3-9) 

n=l 

M 

S.T.  E w = W (3-10) 

i n 

n=l 

wn  > 0 n = 1,  2,  . . M (3-11) 

Note  that  the  objective  function  is  the  sum  of  strictly  convex  functions 
(as  shown  before).  In  addition,  the  feasible  region  is  a convex  set. 

To  solve  this  problem,  we  can  either  eliminate  equation  (3-10)  by 
solving  for  any  of  the  wn's,  or  use  the  Lagrange  multiplier  method. 

The  optimal  solution  for  (RT1),  using  either  method,  is 

^ n 

w^  = — for  n = 1,  2,  . . . , M (3-12) 

and  g*  = FN  (^,  DQI))  f or  n = 1 , 2 , . . . , M (3-13) 

What  does  this  imply  regarding  problem  (Tl)?  If  this  solution  of 
(RT1)  is  feasible  to  (Tl)  [i.e.,  it  satisfies  eqs.  (3-3)  and  (3-4)], 
then  it  is  the  unique  optimal  solution  for  (Tl)  because  the  objective 
function  is  strictly  convex.  Otherwise  there  is  no  feasible  solution  to 
(Tl).  To  prove  the  last  statement,  look  at  eqs.  (3-3)  and  (3-4).  By 
assumption,  at  least  one  of  them  is  not  satisfied,  i.e.,  we  will  have 


FN  (£’  °0D)  + § > ~ EW 


or 
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W 

Since  it  is  not  feasible  to  have  all  wn's  equal  to  — , then  any  other 
feasible  solution  must  have 

w 

w.  > TT  for  some  i 
W 

But  having  w^  > — means  having  a larger  g^  than  before,  i.e., 

Wi  > M * Si  = FN  ^Wi,EW  > FN  ^M,EW 

Thus , we  would  have 

FN  (wi’DoD}  + Wi  > FN  (I»d0D)  + TL  > MT0W  " EW 


or 


W 


FN  (wi »Dq0)  > FN  (h,Dod)  > F 


max 


which  means  no  other  feasible  solution  to  (Tl)  exists, 


In  conclusion,  if  it  is  feasible  to  use  M A/C,  then 


* W 

w =77  for  n = 1 , 2 , 
n M 


• • , M and  g*  - FN  (|, D^) 


Theorem  3-2 


Let  = the  minimum  feasible  number  of  A/C  that  can  be  used. 


a'  D 

If  —p—  > — , then  the  minimum  feasible  number  of  A/C  results 


"2a, 


- M 


min 


in  the  lowest  total  fuel  cost. 

Proof  of  theorem  3-2.  Here  we  need  to  show  that  the  cost  is  an 
increasing  function  of  the  number  of  A/C  used.  First,  note  that  if  M 
A/C  are  used,  then 
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IW 

M n = 1 , 1 , . . . , M 

0 n » M + 1,  . . .,  N 


from  theorem  3-1. 

Let  TF(M)  = the  total  fuel  cost  when  M A/C  are  used  (M _>  1 and  an 
integer).  Thus,  TF(M)  is  a discrete  function. 


M 

TF(M)  = Z gn 
n=l 


M 

L 

n=l 


FN 


n 


(Vdod>  ' 


M * FN 


^ M * °0D  ^ 


+ ^ ^<aif+a;,2+2VoD  (3-i3) 

Let  TF1(M)  = TF(M)  where  M _>  1 and  continuous.  Here,  TF1(M)  is  the 
continuous  version  of  TF(M).  That  is,  it  includes  the  value  of  the 
function  for  integer  as  well  as  noninteger  values  of  M.  Also  both 
functions  should  have  the  same  shape,  (i.e.,  if  TF1(M)  is  increasing,  so 
is  TF(M)).  Thus,  we  need  only  to  show  that  TF1(M)  is  an  increasing 
function  of  M,  or  equivalently  that  the  first  derivative  of  TF1(M)  is 
nondecreasing. 


dTFl(M) 

dM 


-a' 


-rMi  - 


(aiS+a;>  + 


2a, 


OD 


^ ~i  l + a;>2  + 2aiDoD 


(3-14) 


Note  that  a'  >0  and  a,  <0. 
o 1 

Now  mSSl  >0  if  — > 2(a'  ^ + D ) 

dM  - -a1  - o M OD 


a ' D 

. , ,,  O OD  . W 

or  equivalently  if  —x ?—  > — = w 

-2a.  a'  - M n 

1 o 


(3-15) 


Certainly  if  this  is  true  for  M = Mmin,  it  would  be  true  for  M _> 
Minin’  Q.E.D. 
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Let's  see  the  implication  of  this  for  the  C-5A  and  a real-world 
problem.  The  left  hand  side  of  inequality  (3-15)  is  minimum  when  Dqq  is 
maximum  and  still  s-convex,  i.e.,  max  DQD  = 6,220  nautical  miles  (or 
tt/4  in  Radian). 

For  the  C-5A  we  have 

a1  = -0.027,  aQ  = 36.283,  EW  = 320  k lbs 
so  a^  = aQ  + al  EW  = 27.643 


Thus , we  need  to  have 
286  k lbs  > 77  = w 

But  from  the  physical  restriction  on  the  C-5A  we  know  that  wn  <_  200  k 
lbs.  Thus  the  above  condition  is  always  satisfied  for  the  C-5A.  More- 
over is  the  smallest  positive  integer  M that  satisfies  the 

following  two  conditions: 

FN  (^,D__)  < Min  {F  , MTOW  - EW  - ■£} 

M OD  - 1 max’  MJ 

and 


M - 


w. 


max 


w where 
max 

maximum  allowable  cargo  weight  due  to  structural  weight 
limitation 


Figure  3-1  shows  the  total  fuel  cost,  TF(M) , as  a function  of  the 


number  of  A/C  used 
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TF(M)  =Total  Fuel 

in  1000  lbs. 


W = 200,  Doo=5000  mi 


M = Number  of  A/C  used 


Figure  3-1:  Total  fuel  cost  versus  number  of  A/C  used. 


CHAPTER  IV 


PROBLEM  (P2) — MODEL  AND  SOLUTION 


1 A/C 

1 Refueling 
required 

(2  legs) 

Problem  Description 

Consider  the  situation  where  there  is  one  transport  A/C  with  cargo 

weight  wq.  The  fuel  required  to  travel  a distance  DQD — the  distance 

between  the  origin  and  destination — is  FN  (w0,DoD^*  But,  the  maximum 

fuel  the  transport  A/ C can  carry  is  g„„^  (w  ) = Min  fF  , MTOW  - EW  - w 1. 

r J &max  o'  l max  oJ 

So,  if  FN  (w  jD-.p.)  > g (w  ) or  equivalently  if  R[g  (w  ) , w ] < D^p., 

’ o’  OD'  6max  o'  M J 6max  o’  oJ  OD* 

then  the  transport  A/ C must  refuel  en  route  in  order  to  reach  its  desti- 
nation. Now,  assume  there  is  a third  base  that  has  a tanker  A/C  ready 
to  refuel  the  transport  A/C  in  mid-air.  The  maximum  fuel  the  tanker  can 
carry  is  H x#  In  this  case  both  A/ C would  fly  to  a point  (g,^)  where 

they  meet.  There,  some  fuel  would  be  transferred  from  the  tanker  to  the 

transport  A/C.  The  tanker  then  would  return  to  its  base,  and  the  trans- 
port A/C  would  continue  to  its  destination.  The  objective  is  to  deter- 
mine the  location  of  the  refueling  point  so  that  the  total  fuel  require- 

ment for  all  A/C  is  minimized.  There  are  some  limitations  on  the 
location  of  the  refueling  point.  To  see  this,  let 

d^ (0  ,q> ) = the  distance  between  the  origin  base  (Oq^q)  and 

the  refueling  point  (g  ) 
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d2(M) 

= the  distance  between  the  refueling  point  ( 0,  <j>)  and 
the  destination 

d3(  8, 4>) 

= the  distance  between  the  tanker  base  ( 0 , <j>  ) and 

B B 

the  refueling  point  (0,<t>) 

R (g0’wo) 

= the  range  of  the  transport  A/C  when  its  initial 
fuel  is  gQ  and  cargo  weight  is  wQ 

gmax^wo^ 

= the  maximum  fuel  the  transport  A/ C can  carry  when 
its  cargo  weight  is  wQ  and  still  be  able  to  take- 
off 

GWmax 

= min  {F  , MTOW  - EW  - w } 

max  o 

= the  maximum  allowed  gross  weight  of  a transport 
A/C 

Gmax(wo) 

= the  maximum  fuel  the  transport  can  carry  when  it 

R*(h0,w) 

is  in  the  air,  given  that  its  cargo  weight  is  wQ 

= min  {F  , GW  - EW  - w } 

max  max  o 

= the  range  of  the  tanker  A/C  when  its  initial  fuel 
is  hQ  and  the  cargo  weight  is  w (here,  w = 0) 

FN*(w,d) 

= exact  amount  of  fuel  needed  by  the  tanker  to  fly  a 
distance  d when  its  cargo  weight  is  w (usually 
w = 0) 

FC*(h0,w,d) 

= the  fuel  consumed  by  the  tanker  in  flying  a 

distance  d when  its  initial  fuel  is  hQ  and  cargo 
weight  is  w (usually  w = 0) 

fl 

- ^(M) 

= the  fuel  consumed  by  the  transport  A/C  in  flying 
from  Ce0»<t>0)  to  (0  ,<t>) 

= FC(go,wo,d1) 
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f2  = f2(0»+) 

= the  fuel  needed  by  the  transport  A/C  to  go  from 
(0,<f>)  to  (9D»<j>D) 

= FN(wo,d2) 

^3  = f 3 ( ® » d>) 

= the  fuel  consumed  by  the  tanker  in  going  from 

( 9b«  V to 

= FC*  (ho,0,d3) 

f4  = f4(9»<t>) 

= the  fuel  needed  by  the  tanker  to  go  from  (0,4>) 
back  to  (0g»4>B) 

= FN*  (0,d3) 


There  are  several  conditions  that  must  be  satisfied  by  the  location 
of  the  refueling  point.  First,  the  distance  between  the  origin  base  and 
the  refueling  point  should  not  exceed  the  maximum  range  of  the  transport 


A/C,  i.e.. 


dj^C0  ,<j>)  < RUmax(wo),wo]  (4-1) 

Second,  refueling  has  to  be  done  such  that  the  transport  A/C  is  able  to 
reach  its  destination  from  the  refueling  point,  i.e.. 


d (0 ,4>)  < R[G  (w  ) ,w  ] (4-2) 

2 - max  o o 

Third,  the  tanker  should  be  able  to  make  a round  trip  safely,  i.e., 

d3(e,+)  <}r*(hui,o)  (4-3) 

Also,  there  are  other  conditions  that  must  be  satisfied  by  gQ  and  hQ, 
the  initial  fuel  for  the  transport  and  the  tanker,  respectively.  In 
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fact,  whether  or  not  gQ  and  hQ  are  feasible  depends  upon  the  location  of 
the  refueling  point.  This  is  because  gQ  must  be  at  least  equal  to  the 
amount  of  fuel  needed  to  fly  from  the  origin  to  the  refueling  point, 
i.e. , 


Similarly,  hQ  should  permit  the  tanker  to  make  a trip  from  its  base  to 
( 9,  4>)  and  back,  i.e.. 


Finally,  at  the  refueling  point,  the  amount  of  fuel  left  in  the  tanker 
is  (hQ  - f^)»  and  in  the  transport  is  (gQ  - f^).  To  fly  from  the 
refueling  point  to  the  destination,  the  transport  needs  an  amount  of 
fuel  equal  to  f2«  Thus,  the  amount  of  fuel  to  be  transferred  from  the 
tanker  to  the  transport  A/C  is  = f2  - (gQ  - f^)  = + ^2  ~ ®o*  ^he 

amount  of  fuel  left  in  the  tanker  after  transferring  = (hQ  - f^)  - 
amount  transferred  = (hQ  - f3)  - ( f 1 + f2  “ gQ)  = hQ  + gQ  - fj  - f2  - 
f3.  This  amount  has  to  be  greater  than  or  equal  to  the  amount  of  fuel 
needed  by  the  tanker  to  fly  back  to  its  base,  or: 


FN(wo’dl)  1 ( 9»  <f>)  1 8, 


'max 


(4-4) 


FN*  (0,2d3)  < hQ  (0,*)  < Hmax 


(4-5) 


ho  + So  ~ fl  " f2  " f3  2 f4’  or 


4 


(4-6) 


Formulation  of  (P2) 


The  problem  can  be  stated  mathematically  as  follows:  find  (9,^), 


gQ  and  hQ  that  will 


4 

(P2)  Minimize  V (9  ,<f>  ) = E f ^(9  ,4>  ) 

i=l 


45 


di(  0,  <J>)  < R[gm  (w  ),w  ] 
1 — max  o o 

(4-1) 

d2  (0,4,)  < R[Gmax(wo),wo] 

(4-2) 

d3(e,*)  <^*(Hmax>0) 

(4-3) 

™<Vdl>  < §o  < Smax<"o> 

(4-4) 

F»*(0,2d3>  < h0  < V* 

(4-5) 

4 

^ho  + S0)>  1 f .( 0 f<j>) 

i=l 

(4-6) 

0 _<  0 < 2^ 

(4-7) 

0 4>  £ ir 

(4-8) 

A few  remarks  are  in  order.  First,  the  intersection  of  the  first 
three  inequality  constraints  ((4-1),  (4-2)  and  (4-3))  represent  a region 
inside  which  all  feasible  refueling  points  must  lie.  Figure  4-1  shows 
an  example  of  such  a region.  Note  that  even  though  all  the  feasible 
points  are  included  in  this  region,  not  all  the  points  in  this  region 
are  necessarily  feasible.  For  any  point  in  this  region  to  be  feasible, 
there  must  exist  values  for  its  associated  gQ  and  hQ  that  satisfy 
inequality  (4-6)  which  can  be  rewritten  as 

[ gQ(0  »d> ) ~ f ><t> ) 1 + [ hQ(e  ,<t> ) - f 3(e  ,<{))]  > f2(e,<t>)  + f4(0,<j>)  (4-9) 

The  LHS  of  (4-9)  is  largest  when  both  gQ  and  hQ  attain  their  maximum 
values.  From  this,  it  can  be  concluded  that  any  feasible  point  must 
satisfy  (4-1),  (4-2),  (4-3)  and  the  following  constraint: 
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O Refueling  Point 

Transport  A/C  route 

— — — Tanker  route 


Figure  4-1:  The  region  inside  which  all  feasible  refueling 

points  lie. 


47 


!<W(V  - FCl<WV'Vdl>l  + “Lax  * FC*<Vx-°’d3)1 


>.  f2(0  »4>)  + f4(e ,<)>) 


(4-10) 


which  is  derived  from  (4-4),  (4-5),  and  (4-6).  Therefore,  it  is  easy  to 
check  for  the  feasibility  of  any  point  inside  this  region.  Also,  note 
that  the  RHS  of  (4-6)  is,  in  fact,  the  objective  function.  Therefore, 
if  we  start  with  a feasible  point  and  find  another  one,  inside  the 
region  described  above,  that  has  less  cost,  the  second  point  is  auto- 
matically a feasible  point. 

The  second  remark  also  is  concerned  with  constraint  (4-6).  Since 
values  of  gQ  and  hQ  greater  than  needed  will  increase  the  fuel 
consumption  unnecessarily,  we  can  restrict  our  attention  to  the  case 
where 

4 

g.(0  ,<l> ) + h (6  ,$ ) = 1 f (0  ,4>)  (4-11) 

i=l 

This  means  that  there  is  only  one  value  of  hQ(0,<t>)  for  every  choice  of 
gQ C®  »<J> ) • This  value  can  be  calculated  from  (4-11)  as 

ho(0,4>)  = FN*[hf(gQ),d3]  + hf(gQ)  (4-12) 

where 

hf(g0)  = f 2C©  ,4> ) + f4(0,<l>)  - [gQ(0  »<j> ) - ^(0,*)]  (4-13) 

This  means  that  we  have  to  worry  only  about  selecting  for  each  point  the 
best  value  of  gQ(0  ,<j>  ) (denoted  by  gQ(0,<J>))  that  will  result  in  the 
lowest  total  fuel  cost  for  that  point.  Consequently,  gQ  (0  ,<j> ) is  the 
solution  to  the  following  NLP  subproblem  in  one  dimension  (gQ): 


(SP1)  Minimize  gQ  + FN* [hf (gQ) ,d3 ] + hf(gQ) 


(4-14) 
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S.T. 


FN  (w  , di ) < g < g (w  ) 
o’  1 7 — °o  _ °maxv'wo'' 


(4-4) 


FN*[hf(g0),d3]  + hf(gQ)  £ 


(4-15) 


where  hf(g0)  is  defined  by  (4  13).  This  problem  is  convex  in  gQ  and  can 
be  easily  solved  using  Golden  Section  Search  or  Fibonacci  Search  if  we 
combine  (4-2)  and  (4-15)  into  the  following  single  constraint: 

gmin(9’^  ^ go(0’4>)  < gmax(w0>’  (4-16) 


where 


g«in(  »•  »> 


FN  (w  ,d,) 
o 1 

FN  (gf  + WQ,d1)  + gf 


if  gf(9,  <t>)  < 0, 

(4-17) 

if  gf(  9,  4>)  > 0, 


and 

gf(  0.  <p  - f2(  6,  <t)  + f4(  9,  +)  - [Hmax  - FC*(Hmax,0,d3)]  (4-18) 

The  third  remark  is  that  the  great  circle  distance  does  not  satisfy 
the  usual  convexity  definition  in  En.  However,  this  definition  has  been 
modified  and  extended  to  accommodate  the  situation  where  distance  and 
set  are  defined  on  the  surface  of  a sphere  [1,3,10,12].  The  following 
definitions,  which  describe  spherical  convexity,  are  needed  to  follow 
the  discussion  relative  to  the  convexity  of  the  objective  function  and 
the  feasible  region.  (The  interested  reader  is  referred  to  Appendix  C 
for  a complete  list  of  properties  and  results  of  spherical  convexity.) 


Definition  1 

A set  of  points  on  a sphere  is  said  to  be  spherically  convex  (s- 
convex  set)  if,  for  any  two  points  of  the  set,  the  whole  shorter  great 
circle  arc  connecting  them  is  included  in  the  set  [12]. 
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Note  that  on  a sphere  the  largest  convex  set,  other  than  a whole 
sphere,  is  a hemisphere. 

Definition  2 

A spherically  convex  combination  of  two  points  r^  and  ^ is  defined 
as  the  point  r = p(r^,r^,X)  that  lies  on  the  shorter  great  circle  arc 
between  r^  and  , such  that  the  great  circle  distance  between  r^  and  r 
is  X * d(r  ,r  ) for  0 < X < 1 [12]. 

Definition  3 

A function  f(r)  over  an  s-convex  set  is  said  to  be  a spherically 
convex  function  (s-convex  function)  if  for  every  0 < X < 1: 

fCpCr^.X))  < XfUj)  + (1  - X)  f(r2) 

for  any  element  r^  and  ^ of  the  set  [12]. 

Another  important  result  due  to  Drezner  and  Weslowsky  [12]  says 
that  the  distance  from  a given  point  r is  a s-convex  function  within  a 
circle  of  radius  tt/2  and  center  r.  Since  the  right  hand  side  of 
constraints  (4-1),  (4-2)  and  (4-3)  are  usually  much  less  that  tt/2, 
each  of  these  constraints  represent  a s-convex  set.  The  region 
resulting  from  their  intersection  will  also  be  an  s-convex  set. 

Therefore  Problem  (P2),  without  constraint  (4-10),  will  have  a s-convex 
feasible  region. 

Our  fourth  remark  is  related  to  Aly  et  al.  [1]  who  showed  that  the 
search  for  an  optimal  solution  to  the  spherical  Weber  problem  can  be 
restricted  to  the  spherical  convex  hull  of  the  demand  points.  We  have 
extended  this  result  to  the  case  at  hand  (see  Appendix  D)  to  show  that 
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the  search  for  the  optimal  refueling  point  can  be  restricted  to  the 
spherically  convex  hull  of  the  three  base  points  (the  origin,  the 
destination,  and  the  tanker  base). 

Taking  the  first  and  fourth  remarks  together,  we  can  see  that  the 
search  for  the  optimal  refueling  point  can  be  restricted  to  the  inter- 
section of  the  region  described  by  (4-1),  (4-2),  and  (4-3)  with  the 
spherically  convex  hull  of  the  three  bases. 

Finally,  the  fifth  remark  concerns  the  objective  function.  It  is  a 
nonlinear  function  of  the  spherical  distance  (the  great  circle 
distance).  Using  the  fuel  function  derived  earlier,  it  is  seen  that 

o)2  - 2Vl 


‘o)2  ~ 2ald3 
1 


Vca+a,g0)2  - Zajdj  -Va2+2aj<i2 

’■i(a;+alho>2  - 2ald3  -Vta;)2  + Zajdj 

The  way  the  objective  function  is  written,  it  is  difficult  to 
verify  its  s— convexity  either  by  using  the  second  derivative  test  or  by 


fl  go  + a. 


f~  = - — + 
2 a. 


's/a  + 2a^( 


f~  = h + — 
3 a. 


'Ik  + v 


-a1  yj [a'J"  + 2a,d 


f,  = — 2-  + —2 

4 aL  ax 


13 


I = h + g — 

i=l  1 0 ° ai 


51 


applying  the  definition  of  a s-convex  function  directly.  Nevertheless, 

s-convexity  of  the  objective  function  can  be  shown  by  rewriting  the 

objective  function  using  some  s-convex  function  like  FN(*  ,•  ) . To  do 

this,  define  for  any  point  (6  ^ ) the  following: 

WFL[gQ  (0  ) , d ^ ] = weight  of  the  fuel  left  in  the  transport 

A/ C when  it  flies  a distance  d^  (from  the 
origin  to  (0  ,<)>))  given  that  its  initial 
fuel  is  gQ  (0  ,<{)  ) 

Thus,  WFL  [gQ(0  = gQ(0  ,$)  - FC[go(0  ,4,  ) ,wo,d1]  (4-19) 

where  FN(wQ,d1)  < go(0  ,<> ) < gmax(wQ).  (4-4) 

WFL(*  ,•  ) is  s-convex  because  FC(*  ,•  ,•  ) is  s-concave. 

Let  WFLmax(dL)  = Max  WFL  [gQ(0 

go 

= OTL  [W(wo>*dl]*  (4-20) 

Note  that 

0 < WFL  [ gQ (0  ,<)>  ) , d L ] < WFLmax(d1),  (4-21) 

and  that  WFLfflax(.)  is  s-convex  also. 

•k 

Since  gQ(0  ,<J>  ) is  the  best  value  of  gQ(0,<j)),  then  using  the  s-convex 
functions  FN  and  WFL,  the  elements  of  the  objective  function  can  be 
rewritten  as 


f ^ (0  »4> ) 

= FN  {wo  + WFL[g 

q (0  »<j)  ) >d^  ] , d 

(4-22) 

f 2(9  >4> ) 

= FN(wo,d2) 

(4-23) 

£4(0  >4> ) 

= FN*  (0,d3) 

(4-24) 

f 3 (9  >4> ) 

= FN*  { f 2 + f 4 - 

WFL[g*(0  ,4),d1],d3} 

(4-25) 

4 

and  the  total  fuel  consumption  = V (0  ,<j>  ) = £ f . (0  ,<j>  ) • 

i=l  1 
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Both  FN  (•,•)  and  FN  (•,•)  are  convex  and  increasing  functions  of 
their  arguments,  which,  in  turn,  are  s-convex  functions.  Therefore,  the 
resulting  objective  function  is  s-convex  [17,26].  In  Appendix  D,  we 
present  another  method  to  prove  s-convexity  of  the  objective  function. 
There,  we  give  a sufficient  condition  for  the  s-convexity  which  depends 
on  the  distance  between  the  points  of  the  feasible  region  and  the 
destination.  This  condition  is  arrived  at  using  the  worst  case  method. 


Motivation  for  the  Solution  Procedure  of  (P2) 

For  this  problem,  the  optimal  point  (0*,<j>*)  and  its  optimal  g*(0*,4>*) 

(and  consequently,  its  optimal  hQ(0  ,<J>  ))  are  to  be  found.  Since  gQ( 0, 4>) 

depends  upon  ( 0 , 4>)  and  we  do  not  know  this  in  advance  (remember,  we  have 

to  solve  (SP1)  to  find  that  out!),  it  is  difficult  to  use  (4-22),  (4-23), 

(4-24)  and  (4—25)  directly.  That  is,  we  cannot  search  for  both  (0  , $ ) 

* * * 

and  gQ  ( 9 , <t>  ) simultaneously.  However,  we  can  overcome  this  by  creating 
an  upper  bound  function  that  can  be  improved  from  one  iteration  to 
another  until  its  minimum  value  coincides  with  the  minimum  of  the 
original  function.  To  see  this  let 
WFL  [g_(  0,  <{.)  ,d  ] 

a (9,4>)  wfl  TdT)  (4-26) 

max  l' 

then, 

0 < a ( 0,  <j>)  < 1 (4-27) 


for  all  (0,<j>). 

jSf  ^ ^ 

Thus,  gQ ( 9,  <|>)  would  result  in  a ( 0,  <j>) , where  a ( 0,  <j>)  is  the  best  value 
of  a ( 0,  <j>) , i.e. , 


* WFL[g  (0,4>),d 

a (0,4>)  « ° 1 


WFL  (dT) 
max  1 


(4-28) 
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Also,  from  (3-38),  we  have 


WFL  [g  (9,<(>)  ,d  ] = oCM)  * WFL  (d.)  (4-29) 

° 1 max  1 

The  element  of  the  objective  function  can  be  rewritten  using 
[a  (9,<j>)  • WFLmax(d^]  instead  of  WFL[gQ(  9 , <)>)  ,d^  ] as: 

f 1 (e  ,4>)  = FN  [w  + o*(9,4>)  • WFL  (d.  ) ,d  ] 

f2(M)  = FN  (wo,d2) 

f4(0 ,4>)  = fn*  (0,d3) 

f3(0,4>)  = FN*  [f2  + f4  - ct* ( 9 » 4> ) • WFLrflax(d1),d3] 

Now  setting  a(0,<)>)  = a for  all  points,  and  denoting  the  resulting 
function  by  UB(9,<)>,a),  one  gets 


-©- 

CD 

•— 1 

o 

s 

S3 

II 

WFWdi>-dil 

(4-30) 

f 2 (0  »4>) 

= f 2 ( 0 » d5 ) 

- F»  C»0,d2> 

(4-31) 

»♦) 

= f4(e,4>) 

- FN*  (0,d3) 

(4-32) 

f 3 (0  »4>) 

= fn*  {f2 

+ £ - a • WFL  (d  ),d  } 

4 max  1 3 J 

(4-33) 

where 

0 < a < 

1 

and  UB(9  ,<j>  ,a) 

4 

= E 

i=l 

fJ.(M) 

(4-34) 

This  function  has  the  following  characteristics: 

1)  It  is  s-convex  because  it  is  the  sum  of  s-convex  functions; 

2)  UB[9  ,4>  ,o ] > V(0  ,4>)  (4-35) 

In  other  words,,  since  any  g^  other  than  g (0,<j>)  would  result  in 

o o T 

higher  cost  for  V(0  ,4> ) , and  also  produce  an  a that  is  different 
* 

from  a (0  ,<J> ) , then  any  a used  in  UB(9,<|>,a)  would  make 

UB(9  ,<J>  ,a)  > V(0  ,4> ) . Therefore  UB  is  an  upper-bound  for  the  objective 

function; 


54 


(4-36) 


*Jg 

That  is  the  equal  sign  of  eq.  (4-35)  holds  true  when  a = a (e,^). 
Therefore 


UB[9  ,<J>  ,a  (0  ,<j>  )]  = V(0  , ) 


* * * * * * * 


(4-37) 


4)  UB(0,<j),a)  = V ( 0 , (j>)  at  the  boundary  points  where  d (9,^)  = 

R[gmax(wo)  ,Wo^  because , at  those  points,  WFL]nax(d1)  = 0.  Thus  the 
value  of  a does  not  matter. 

How  can  we  use  this  UB  function  to  solve  (P2)?  Given  a 
particular  value  of  a,  we  can  optimize  the  UB  function  to  find  the 
best  ( 0 , <j>)  that  corresponds  to  it.  Now,  if  at  some  iteration  we 
have  the  value  of  a (9  ,<(,  ) [without  knowing  (9*,^,*)],  and  we 

jlf  jjg 

optimize  the  UB  function,  then  we  will  end  up  with  (9  , ^ ) . This  is 
true  since  if  (0  ,,)>  ) minimizes  the  true  objective  function,  then  it 

jflg  "ft 

will  also  minimize  UB  [ 9 , (j,  >ct  (9  , ^ )]. 


Solution  Procedure  of  (P2)  and  Convergence  Proof 


1 .  Initialization: 

Let  k = 0,  and  a°  = 1 . Start  with  a point  (0O,<j,°)  that  satisfies 
(4-1),  (4-2)  and  (4-3)  (Selection  of  a good  starting  point  is 
discussed  on  page  56). 


2.  Let  k - k + 1.  Find  the  point  that  will  minimize  UB( 9 , (j,  ,ct^C  S s.t. 
(4-1),  (4-2),  and  (4-3).  Let  the  solution  be  (9k,<j)k). 


3.  For  (9  ,^R) , solve  (SP1)  to  find  the  best  g^g^,^).  Use  equation 
(4-28)  to  find  ci*(0k,<t>k). 


* k k 
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* k k k- 1 

4.  Ifa(9,<j>)*a  ,we  have  an  improvement  in  the  objective  func- 
tion. Therefore,  let  ak  = a (0k,<(>k)  and  go  to  step  2.  Otherwise, 
the  solution  is  (0  ,<j>  ),  gQ(9k,<j)k),  and  1^(0*, <j>K)  is  calculated  from 
(4-12)  and  (4-13).  Stop. 


Next,  a convergence  proof  of  this  procedure  is  given  below. 


We  need  to  show  that 


£ fi(ek“1,<|>k“1)  > z f (9k,<j>k) 

i=l  i=l  1 


(4-38) 


Assume,  at  iteration  k-1 , that  a*(0k-1  ,(j>k_1 ) * ak~2.  (Otherwise, 
we  would  have  stopped  there.)  By  setting  ak-i  = a*(0k_1  ,<j)k-1 ) , we  have 


Z f (0k-1,$k-1)  = UB(0k-1 ,4>k— 1 ,ak_1 ) 
i=l 

At  iteration  k,  since  (0k ,4>k)  minimizes  UB(0  ,<j>  ,ak_1 ) , then 
UB(0k"1,4,k"1,ak_1)  > UB(0 k ,<j>k ,ak_  1 ) 

But , 

UB(9k,(j>k,ak_1)  > Z f . (9k,4>k) 

“ i=l  1 

and  inequality  (4-38)  follows. 

Moreover , 


(4-39) 


(4-40) 


(4-51) 


4 

5 

i=l 


(9k_1,4>k_1; 


4 

> Z 

i=l 


fi(9k,<|)k). 


if  UB(0k,<j>k,ak-1 ) > 


4 

Z 

i=l 


(4-52) 


(3-55) 
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The  last  inequality  holds  true  if  a*(0k,<]>k)  * ak_1.  Therefore,  the  true 
objective  function  improves  from  one  iteration  to  another.  Since  the 
problem  is  s-convex,  the  solution  procedure  will  converge  to  the  global 
optimum  solution  [12].  0 E D 


Finding  a Good  Starting  Point 

Here,  the  following  lower  bound  function  is  used  to  find  a good 
starting  point. 

= FN(w  d ) + FN(w  ,d„)  + FN* [h  , (0,<j,),dj 
o 1 o L f,min  3 

+ FN*(0,d3)  (4-44) 

where 

hf,min(0,<f,)  = FN(w0’d2)  + FN*(°»d3) 

“ {g  (w  ) - FC  [g  (w  ) ,w  ,d,  ] } 

1 max  o max  o o 1 > 

= FN(wQ,d2)  + FN*(0,d3)  - WFLmax(d1 ) (4-45) 

This  lower  bound  function  is  s-convex,  and  LB(Q  ,4, ) < V(0,<j>).  A physical 
explanation  for  LB(0,<}>)  < V(0,<j>)  follows: 

Each  term  in  the  LB  function  represents  the  minimum  amount  of  fuel  the 
various  A/C  need  for  each  particular  leg  regardless  of  the  actual  fuel 
consumed  or  needed  for  the  other  legs.  That  is,  for  each  leg  of  the 
flight  we  assume  the  other  A/C  will  come  to  the  refueling  point  with  the 
maximum  amount  of  fuel  possible.  The  equality  holds  at  the  boundary 
points  where  d (0,<j>)  = R[gm  (w  ) ,w  ].  How  can  the  LB  function  be 
used?  Start  with  any  point  inside  the  region  described  by  the 
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intersection  of  inequalities  (4-1),  (4-2)  and  (4-3),  and  solve  the 
following  problem: 

Min  LB(q  ,(j, ) 

S.T.  (4-1),  (4-2)  and  (4-3) 

If  the  optimal  solution  to  LB(0  ,4)  lies  on  the  boundary  where  d^  = 
R[gmax(wo’wo]’  then  iC  is  the  optimal  solution  to  (P2),  and  we  do  not 
have  to  search  any  further.  If  the  minimum  value  of  the  lower  bound 
function  is  > Hfflax  + gmax(wQ),  then  there  is  no  feasible  solution  to 
this  problem. 

There  is  no  guarantee  that  the  resulting  point  will  be  feasible  if 

Min  { LB(0  ,4,)}  < H + g (w  ); 

(9,4)  max  max  0 

however,  its  feasibility  could  be  checked  using  inequality  (4-10).  If 
it  is  feasible,  then  all  subsequent  points  will  be  feasible  also.  (See 
the  first  remark  about  the  objective  function  on  page  45.)  Even  if  it 
is  not  feasible,  this  will  not  render  the  solution  procedure  of  (P2) 
useless  because  the  feasibility  of  the  starting  point  is  not  necessary 
as  long  as  the  point  satisfies  inequalities  (4-1),  (4-2)  and  (4-3).  The 
procedure  would  still  be  valid  because  the  objective  function  improves 
from  one  iteration  to  another.  Improvement  and  convergence  do  not 
depend  on  the  point  satisfying  inequality  (4-10).  In  fact,  we  need  to 
check  the  feasibility  of  the  point  only  at  the  end  when  the  procedure 
stops  and  convergence  is  reached.  We  then  check  to  see  if  the  final 
point  satisfies  inequality  (4—10).  If  it  does,  then  the  problem  has  a 
feasible  solution  and  the  solution  at  hand  is  the  optimum  feasible 
solution  we  are  seeking.  Otherwise  the  problem  does  not  have  a feasible 


solution  at  all. 


58 


Computational  Results 

The  above  mentioned  procedure  was  programmed  on  a Vax  11/750  using 
BASIC.  There  was  no  attempt  to  make  it  efficient.  Several  test 
problems  were  solved  to  represent  a variety  of  possible  geographic 
configuration  and  complexities.  The  minimum  CPU  time  was  approximately 
0.5  second  (2  iterations)  and  the  maximum  was  almost  1.0  second  (4 
iterations)  for  convergence  within  one  degree  longitude  and  latitude. 
Table  4-2  shows  these  test  problems  with  their  results,  while  Table  4-1 
provides  the  key  to  understand  the  entries  of  Table  4-2. 
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TABLE  4-1:  Key  for  trial  runs  chart  of  Table  4-2 


Location 

Longitude 

Latitude 

Code 

New  Jersey 

75W 

40N 

A 

Delaware 

75W 

38N 

B 

North  Carolina 

78W 

35N 

C 

Puerto  Rico 

66W 

18N 

D 

Azores  Islands 

25W 

37N 

E 

Iceland 

20W 

65N 

F 

Germany 

10E 

50N 

G 

Turkey 

30E 

40N 

H 

Saudi-Arabia 

47E 

25N 

I 

Egypt 

28E 

30N 

J 

England 

OE/W 

52N 

K 

W = Cargo  Weight  Code 

100.000  lbs  1 

200.000  lbs  2 

Fuel  ACCR  = Fuel  Accuracy 

100  means  to  within  100  lbs 
10  means  to  within  10  lbs 

POS  ACCR  = Position  Accuracy 


1 means  to  within  1 degree 

2 means  to  within  1 minute 

3 means  to  within  1 second 


TABLE  4-2:  Trial  runs  for  Problem  (P2) 
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CHAPTER  V 


PROBLEM  (P3)—  MODEL  AND  SOLUTION 
Two  A/C 

1 Refueling  Required 
(2  Legs) 

Formulation  of  (P3) 

Consider  the  case  where  there  are  two  transport  A/C  available  at 
the  origin  base  "0,"  and  two  tanker  A/C  at  another  base  ready  to  refuel 

the  transport  A/ C if  necessary.  It  is  desired  to  find  out  how  to  divide 

the  cargo  weight  optimally  among  both  transport  A/C.  Note  that  the 
chosen  value  of  Wj^  might  necessitate  an  aerial  refueling,  in  which  case 
the  optimal  refueling  point  also  needs  to  be  determined.  Also,  it  may 
not  be  necessary  to  use  both  A/C.  Rather,  we  may  use  only  one  A/C  if 
feasible.  The  chosen  option  will  be  the  one  which  results  in  the  lowest 
total  cost.  To  find  this  we  need  to  solve  the  following  problem: 

(P3)  TC2(W)  = Min  TC(w1,w2)  = Min  [C(wx)  + C(w2)j 
S.T.  + w2  < W 

wj  > 0 , w2  >0 

where,  TC2(W)  = Total  fuel  cost  incurred  by  transporting  cargo  which 
weighs  W when  there  are  only  two  A/C  available,  and  C(w)  = cost  of 
transporting  a cargo  weight  of  w by  one  A/C. 


61 


62 


10  if  w = 0 

C0(w)  = FN  (w,D0D)  if  0 < w < Zq 

C, (w)  = solution  of  of  (P2)  if  Z < w < Z, 

1 o 1 

00  if  ZL  < w 

TC(w-^,W2)  = C(w^)  + C(w2),  where 

Zo  = the  raaximum  cargo  weight  one  transport  A/C  can  carry  from  "0"  to 
"D"  without  refueling. 

Z1  = the  maximum  cargo  weight  one  transport  A/C  can  carry  from  "0"  to 
"D”  by  refueling  once.  (Z1  > ZQ.) 

At  the  end  of  this  chapter  we  present  simple  ways  to  compute  ZQ  and 
Zp  Figure  5-1  shows  a typical  C(w)  function.  This  function  is  piece- 
wise  convex  and  strictly  increasing.  The  jump  at  ZQ  depends  on  DQD,  W 
and  the  location  of  "B." 


Solution  of  (P3) 

To  make  the  solution  of  (P3)  easier,  we  can  rewrite  it  in  terras  of 
one  variable  instead  of  two,  i.e., 

(P3)  TC2(W)  = Min  TC(w1 , W - = Min  CCwj)  + C(W  - 

S.T.  0<w^<W 

Solving  this  problem  allows  us  to  identify  seven  (7)  cases,  depending 
upon  the  value  of  W relative  to  ZQ  and  Z^.  These  cases  have  simple 
solutions.  The  analysis  of  these  seven  cases  is  presented  next, 
followed  by  a summary  of  a simple  solution  procedure. 
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C(w)  using  1 A/C 


w = cargo  weight  in  1000  lbs. 


Figure  5-1:  Cost  when  refueling  is  used 
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Case  1:  (ID,  2D)  W < Z 

. f o 

Here,  one  A/C  can  carry  the  whole  cargo  directly  to  the  destination 
without  refueling.  The  cost  is  CQ(W)  = FN  (W,DQD).  Also,  this  implies 
that  two  A/C  can  carry  the  cargo  without  refueling.  From  problem  (PI), 
theorem  (2-1),  we  know  that  it  is  best  to  divide  W equally  among  both 
A/C.  The  cost  will  be  2 CQ  (— ) = 2 FN(^,  D^).  This  has  to  be  compared 
with  CQ(W),  and  the  option  with  the  least  cost  should  be  chosen. 

Case  2;  (1R,  2D)  ZQ  < W < Min  {2  Zo,Z,} 

Again  two  A/C  can  carry  the  whole  cargo  directly  to  the  destina- 
tion, but  aerial  refueling  is  needed  if  only  one  A/C  is  used.  To  find 
out  how  much  the  one  A/C  scenario  costs,  we  need  to  solve  problem  (P2) 
for  w0  = W and  also  find  the  optimal  refueling  point.  The  cost  is  C^(W) 
= V(Q  ,<(>).  Again  costs  for  both  scenarios  have  to  be  compared.  Figure 
(5-2)  shows  TC  (w^ , W - w^)  as  a function  of  w^  for  this  case. 

Case  3:  (2D)  Z,  < W < 2 Z 

i o 

Here  one  A/ C can  not  transport  the  whole  cargo  (W)  even  if  aerial 

refueling  is  used.  But,  two  A/ C can  transport  the  cargo  directly.  The 
W 

cost  is  2 CQ(y) , and  this  is  the  only  option  we  have. 

Case  4:  (1R,  2R,  1D1R)  2ZQ  < W <;  Z1 

Here,  two  A/C  cannot  carry  the  whole  cargo  (W)  directly  to  "D."  At 
least  one  of  them  has  to  refuel.  If  only  one  A/C  refuels,  then  the 

•fa  A 

optimum  weights  are  = Zq  and  w2  = W - ZQ.  The  cost  is  CQ(Zo)  + C^W 

- Z )•  Therefore,  problem  (P2)  has  to  be  solved  for  w = W - Z to 
u o o 

find  C^(W  - ZQ)  and  the  optimum  refueling  point.  But  if  both  A/C 
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TC  (w-j , W-Wl) 


wn  = cargo  weight  of  first  transport  A/C  in  1000  lbs 


w 


*_W 
1 ” 2 


Figure  5-2:  Cost  for  Case  (2) 
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refuel  then  the  optimal  solution  is  w^  = w2  = W/2,  and  the  total  cost 
is  2 C^(W/2) . Thus,  we  have  to  solve  problem  (P2)  for  wq  = W/2  to  find 
the  optimal  refueling  point  and  the  value  of  C^W/2).  Also,  since 
W<  Z^  , then  one  A/C  would  be  able  to  carry  the  whole  cargo  when  aerial 
refueling  is  used.  Again,  we  need  to  solve  problem  (P2)  for  wQ  = W to 
find  CX(W)  and  the  optimal  refueling  point.  Finally,  all  three  costs 
have  to  be  compared  and  the  least  expensive  option  is  chosen.  Figs.  (5-3) 
and  (5-4)  show  a typical  TC(w^,  W - w^)  curve  as  a function  of  w^  for 
this  case. 

Case  5;  (2R,  1D1R) 

2Z  < W < Z + Z,  and  Z,  < W -*•  Max  { 2Z  , Z,}  < W < Z + Z 
O Oi  x lo  lJ  ol 

Again,  two  A/ C have  to  be  used,  and  at  least  one  of  them  must 
refuel.  The  "two  A/C  refueling"  option  costs  2 C^W/2)  which  is 
obtained  by  solving  problem  (P2)  for  wQ  = W/2.  The  other  option,  where 
one  A/C  refuels  and  the  other  flies  directly,  costs  CQ(Z0)  + C^(W  - 
ZQ).  The  cost  C^(W  - ZQ)  is  found  by  solving  (P2)  for  wQ  = W - ZQ. 

Again,  both  costs  have  to  be  compared. 

Case  6:  (2R)  ZQ  + Z1  < W < 2ZX 

Here,  the  two  A/C  have  to  be  used,  and  both  have  to  refuel.  This 
is  the  only  option.  Thus,  the  cost  is  2 C^(W/2),  and  (P2)  has  to  be 
solved  for  wQ  = W/2  to  find  this  cost. 
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TCfw-j.W-v^) 


w1  = Weight  of  cargo  in  first  transport  A/C 
w* * Z0 


Figure  5-3:  Cost  for  case  (4)  2ZQ  < W < Z^. 
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w1  = Weight  of  cargo  in  first  transport  A/C  in  1000  lbs 


Figure  5-4:  Another  example  for  the  cost  of  case  (4) 
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Case  7:  (infeasible  problem)  2Z1  < W 

Here,  two  A/ C are  not  enough  to  carry  the  whole  cargo  even  If  they 
are  both  fully  loaded  and  mid-air  refueling  is  used.  Therefore,  more 
A/C  are  needed,  or  aerial  refueling  might  be  used  twice. 

Before  we  present  the  summary  of  the  solution  procedure,  we  now 
offer  a simple  way  to  find  ZQ  and  Zy 

Finding  the  Value  of  Z 

_ o 

1.  Solve  the  following  equation  for  w:  MTOW  - EW  - w = FN  (w,  DQD), 

where  FN  (w,  DQI))  is  defined  in  eq.  (2-20).  Call  this  value  w . 

3l 

2.  Solve  the  next  equation  for  w:  F^  = FN  (w,  D0D) . Call  this  value 

wb* 

3.  Set  Zq  = Min  jw^,  w^,  w^^j  where  wmax  = the  maximum  weight  that  can 
be  carried  inside  the  A/C  due  to  structural  weight  limitations. 

This  procedure  is  based  on  the  fact  that  Zq  must  satisfy  the 
following  conditions: 

Zo  * wmax’ 

™ <Zo>  °0D)  “ 8„(7'o)  ■ "here 

fWV  * OT°"  - ™ - zo}. 

Finding  the  Value  of  Z] 

The  real  issue  is  determining  if  wQ  < Z^,  for  any  value  of  wQ. 

That  is  we  need  to  see  if  W < Zj,  W/2  < Zj  or  W - ZQ  < Zj  (or  any  other 
value  of  wQ  for  that  matter).  This  is  accomplished  by  solving  problem 
(P2)  for  that  particular  value  of  wQ  and  determining  if  V(9,<j>)  < Hmay  + 
smax  (V-  this  relation  bolds  true,  then  the  problem  is  feasible 
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and  wQ  < Zj.  As  a by-product  we  get  the  value  of  the  total  fuel 

consumption  C1(wq)  which  is  equal  to  V(  9,  <j»)  , together  with 
* * * * 

h.  atld  (9  , <j>  ).  Note  that  Z.  cannot  exceed  the  value  of  w , thus 
u u * max ’ 

we  do  not  need  to  solve  (P2)  unless  wQ  < w(nax.  In  other  words,  we  can 

say  that  w„  > Z,  if  w > w . 

o 1 o max 

Summary  of  Solution  Procedure  of  (P3) 

1 . Find  Z , 

o 

2.  If  W > ZQ  go  to  step  (3).  Otherwise  compare  CQ(W)  with  2 CQ(W/2). 
If  C0(W)  < 2 CQ  (W/2),  use  one  A/C  flying  directly  (ID).  Its 
iriitial  fuel  is  C^(W).  Otherwise  use  two  A/ C (2D),  each  having 
initial  fuel  equals  to  CQ(W/2).  Stop. 

3.  If  W > 2 ZQ,  go  to  step  5.  Otherwise,  solve  (P2)  for  wQ  = W to  see 

if  W < ZL.  If  V ( 0 , 4>)  > Hmax  + gmax  (w),  go  to  step  4 because  W > Z1 

(That  is  (P2)  with  wq  = W is  infeasible).  Otherwise,  compare  C^W) 

= V(  9 , <|>)  with  2 C0(W/2).  If  CX(W)  < 2 CQ(W/2),  use  one  A/C  (1R). 

Get  the  values  of  gQ,  hQ  and  (9  , 4>  ) from  the  solution  of  (P2). 
Otherwise  use  two  A/C  flying  directly  (2D),  where 

£ -k 

gl  = g2  = CQ(W/2)*  Stop* 

4.  Use  two  A/C  flying  directly  (2D),  g*  = g*  = CQ(W/2).  Stop. 

5.  Solve  (P2)  with  wq  = W/2  to  see  if  j < Z . If  (P2)  is  feasible  go 

to  step  6.  Otherwise,  we  have  W > 2 Z^  which  means  two  A/C  cannot 
carry  the  whole  cargo  even  if  aerial  refueling  is  used.  Therefore, 
we  need  more  A/ C or  need  to  refuel  twice,  and  problem  (P3)  is 
infeasible.  Stop. 

Solve  (P2)  , again,  but  for  wQ  = w - ZQ  to  see  if  w < ZQ  + Z1  (or  W - 
Zo  * ZP*  the  new  (p2)  is  feasible  go  to  step  7.  Otherwise,  we 


6. 
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use  two  A/C,  both  refueling  (2R).  w*  = w*  = W/2  . g*,  h*  and  (0*  <fc*) 

12  o o Y 

are  obtained  from  the  solution  of  the  first  (P2)  problem  mentioned  in 
step  5.  Stop. 

7.  Solve  (P2),  this  time,  for  wq  = W to  see  if  W<  If  the  recent 

(P2)  problem  is  feasible,  go  to  step  8.  Otherwise,  compare  the  cost 

of  step  5 and  step  6,  i.e.,  compare  2 C, (W/2)  with  C (Z  ) + C,  (W2  - 

l o o 1 

Zq).  If  2Cl(W/2)  is  lower,  we  have  = w/2  and  the  two  A/C 

will  need  to  refuel  (2R).  Get  the  values  of  g*,  h*  and  (9 *,<j>*)  for 
both  A/C  from  the  solution  obtained  in  step  5.  Otherwise,  we  have 
= ZQ  an<*  w2  = W ” Z0  Get  the  values  of  g*,  h*  and  (0  *,<(>*) 

for  the  second  A/C  from  the  solution  obtained  in  step  6.  Stop. 

8.  Here,  we  have  the  following  three  options: 

(a)  wx  = W,  w2  = 0,  with  cost  = C^W)  obtained  in  step  7.  (1R) 

(b)  wx  = Zq,  w2  = W - ZQ,  with  cost  = CQ(Z0)  + CL(W  - ZQ)  obtained 

in  step  6.  (1D1R) 

. * * 

Cc)  wx  = w2  = W/2,  with  cost  = 2 C^W/2)  obtained  in  step  5.  (2R). 

Compare  the  cost  and  pick  the  lowest  one.  Get  the  values 
* * * * 

of  gQ»  (0  ,<j>  ) from  the  solution  obtained  in  the  equivalent 
step.  Stop. 

Figure  5-5  gives  a graphical  summary  of  the  solution  procedure, 
while  Fig.  5-6  shows  how  the  fuel  cost  varies  with  W,  the  total  cargo 
weight,  for  a C— 5 A.  Also  it  shows  how  W will  be  transported. 

Before  we  end  this  chapter,  we  would  like  to  give  an  example  to 
illustrate  how  simple  the  aforementioned  solution  procedure  is. 

Example  1: 

Suppose  we  have  2 (identical)  A/C,  and  that  ZQ  = 80  k lbs  and  Z1  = 

180  k lbs  (known).  What  would  the  solution  be  for  different  values  of 
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(2R) 


Figure  5-5:  Graphical  summary  of  the  solution  procedure  of  (P2) 
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Figure  5-6:  Fuel  cost  as  a function  of  W for  a typical  C5-A. 
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W?  Specifically,  what  is  the  solution  of  W = 230  k lbs.?  What  happens 
if  we  do  not  know  the  value  of  Z^? 

Solution: 


I.  Table  5-1  shows  the  solution  when  W lies  in  any  interval. 

II.  For  the  particular  value  of  230  k lbs  for  W,  we  have  two  options 
to  choose  from: 

a)  To  have  2 refueling  A/ C with  wj  - w2  = W/2  = 115  k lbs.  and 

total  fuel  consumption  = TCa  = 20^(115).  The  value  of  0^115) 
is  found  by  solving  (P2)  for  wQ  = 115. 

b)  To  have  one  A/ C fly  directly  to  the  destination  with  w1  = ZQ  = 

80  k lbs.,  and  the  other  to  have  w2  = W - ZQ  = 230  - 80  = 150 
k lbs.  The  fuel  consumption  for  this  option  = TCb  = CQ(80)  + 
0^(150).  The  value  of  0^(150)  is  found  by  solving  (P2)  for  wQ 
= 150.  Finally,  we  compare  TCa  and  TCb  to  choose  the  least 
cost  option. 

III.  Now,  suppose  we  are  not  given  the  value  of  Z ^ in  advance.  We  have 
to  proceed  in  the  following  way  as  suggested  by  Fig.  5-5.  Note 
that  W = 230  > 2Zq  = 160.  Now,  we  have  to  solve  (P2)  for  W/2  = 

115  to  see  if  W/2  < Z^.  When  we  solve  this  problem,  we  see  it  is 
feasible  and  we  also  get  the  value  of  C^W/2  = 115).  Next,  we 

check  if  W - ZQ  = 230  - 80  = 150  < Z1  by  solving  (P2)  for  W - Zq  = 

150.  We  find  out  that  it  is,  and  we  also  get  the  value  of  C^W  - 
ZQ  = 150).  Finally  we  check  if  W = 230  < Zj  by  solving  (P2)  for  W 
= 230.  We  find  that  it  is  not.  This  means  we  have: 

MaxjZ^  2Zq}  < W < Zq  + Z . 
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TABLE  5-1:  Solution  of  example  1 


Interval  for 
W (k  lbs) 

0 

to 

80 

80.001 

to 

160 

160.001 

to 

180 

180.001 

to 

260 

260.001 

to 

360 

above 

360 

ID 

1R 

1R 

2R 

2R 

Need 

or 

or 

or 

or 

more 

Solution(s) 

2D 

2D 

2R 

1D1R 

A/C 

or 

or 

1D1R 

refuel 

twice 
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The  solution  would  be  the  best  of  the  two  following  options 
option  a (2R)  with  cost  = 20^(115) 
and  option  b (1D1R)  with  cost  = 0^(80)  + <^(150) 


as  before. 


CHAPTER  VI 


PROBLEM  (P4) — MODEL  AND  SOLUTION 
N A/C 

1 Refueling  required 
(2  Legs) 

Problem  Description 

Consider  the  situation  where  there  are  N transport  A/C  at  the 
origin  base  0.  Some  or  all  of  these  A/C  can  be  used  to  transport 
cargo  weighing  W to  the  destination  base  "D."  It  is  desired  to  know  how 
many  of  these  A/C  to  use  and  how  to  divide  the  cargo  weight  among  them 
optimally  in  order  to  minimize  total  fuel  consumption.  If  the  load  of 
some  of  these  A/C  exceeds  ZQ,  then  aerial  refueling  is  necessary  and  we 
also  have  to  find  the  optimal  refueling  point. 

Solution  of  (P4) 

The  solution  procedure  is  an  extension  of  the  procedure  used  to 
solve  Problem  (P3).  First,  we  have  to  find  the  different  feasible  ways 
in  which  we  can  transport  the  cargo.  Then,  compare  their  costs  and 
choose  the  one  that  has  the  lowest  total  cost.  To  find  these  combina- 
tions we  first  find  the  values  of  Mm^n  and  Q which  are  defined  as 

n = the  minimum  number  of  transport  A/C  needed  to 

transport  the  cargo  directly  without  aerial 
refueling. 

Q = the  minimum  number  of  transport  A/C  needed  to 

transport  the  cargo  when  aerial  refueling  is  used 
for  all  of  them. 
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Next  we  show  how  we  can  find  Mmin  and  Q.  Also,  we  will  give  an  example 
to  illustrate  the  procedure  of  finding  the  different  feasible  solutions. 
Finally,  a summary  of  the  solution  procedure  is  given. 


Finding  Mm.n 


M 


[-f-1  where 


min 


(6-1) 


f X 1 - the  smallest  integer  greater  than  or  equal  to 


Finding  Q 


q = HH 


(6-2) 


If  is  not  known  then  we  can  find  Q as  follows: 


1.  Since  Z1  < wmax,  let  q = [“4  . 


max 


2.  Solve  problem  (P2)  with  w_  = — to  see  if  it  is  feasible. 


3.  If  it  is  feasible,  then  set  Q = q and  stop.  Otherwise  make  q = q + 

1 and  go  to  step  (2). 

Example 

Suppose  we  have:  W = 490,  ZQ  = 120,  Zj  = 180  and  w^  = 200.  Find 

the  different  feasible  solutions 


w r 490 

^min  I no  1 ^ 


q = r — i = 3 
y 1 180  1 J 


(even  if  we  do  not  know  , we  can  find  Q as  suggested  before,  i.e.. 
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Since  q - ["— - — ] = 3,  solving  problem  (P2)  with  w = — = 163.33 
max  ° 9 

we  find  that  q is  feasible.  Thus  Q = q = 3.) 

So  far  we  have  two  feasible  solutions,  namely  (5D)  and  (3R),  i.e., 

five  A/C  can  carry  the  whole  cargo  directly  or  three  A/C  can  carry  it 

but  they  need  to  refuel  in  mid-air.  Now  let  us  see  if  [ID  2 R]  is 

feasible.  We  get  this  combination  from  (3R)  by  reducing  the  number  of 

refueling  A/C  by  one,  and  letting  one  A/C  fly  directly.  Since  we  can 

have  up  to  ZQ  in  the  A/C  that  is  flying  directly,  we  then  check  to  see 

if  (Q— 1 ) refueling  A/C  can  carry  (W-Zq)  weight.  That  is,  we  check  if 

two  refueling  A/ C can  carry  370  k lbs.  To  do  this  we  see  if 

W - Z 

— < Z 

Q - 1 1 » 

370 

(i.e.,  check  whether  - < 180).  Since  it  is  not,  we  check  whether  [2D 

(Q  - 1)  R]  is  feasible,  i.e.,  check  whether  (2D  2R)  is  feasible.  Again, 
this  translates  to  determining  if 
W - 2Z 


Q - 1 


< zlf 


250 

(i.e.,  whether  - < 180.)  Since  it  is,  we  now  have  a new  feasible 
solution  which  is  (2D  2R).  Starting  with  this  new  solution,  we  try  the 
same  process  of  reducing  the  number  of  refueling  A/C  by  one  and 
increasing  the  number  of  A/ C that  are  flying  directly  by  one  and  check- 
ing its  feasibility.  If  not  feasible,  then  again  increase  the  number  of 
A/C  that  are  flying  directly.  Thus,  we  check  whether  (3D  1R)  is 
feasible  by  determining  if 
W - 3Z 

o 130 

q _ 2 * zi  or  in  our  example  if  — j—  < 180? 
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Since  this  is  true,  we  have  another  new  solution,  namely  (3D  1R).  Now 

we  repeat  this  to  see  if  (4D)  is  feasible.  But  we  know  that  M . =5 

mm 

implies  that  (4D)  is  not  feasible  and  therefore  the  feasible  solutions 

are:  (5D),  (3D  1R) , (2D  2R) , (3R). 

Now  we  come  to  the  problem  of  how  to  divide  the  cargo  among  the  A/C 

to  be  used.  For  the  solution  (5D)  and  (3R)  [or  equivalently  (M^  D) 

and  (QR)],  we  divide  the  cargo  weight  equally  among  the  A/C.  This  is 

because  in  both  cases  the  total  cost  function  is  continuous,  increasing 

and  convex  in  their  respective  region.  This  is  also  shown  to  be  optimal 

in  problems  (PI)  and  (P2).  Now  we  come  to  the  solution  (3D  1R)  and  (2D 

2R) . We  saw  from  the  solution  of  problem  (P3)  that  whenever  we  have  a 

combination  of  A/ C of  which  some  fly  directly  to  the  destination  while 

the  rest  refuel,  the  cargo  load  for  the  A/C  flying  directly  should  be 

* 

set  at  maximum.  That  is,  w^  = for  A/ C flying  directly  while  the  rest 
of  the  cargo  weight  is  divided  equally  among  the  refueling  A/C.  In 
conclusion,  the  different  solutions  for  the  example  problem  are  as 
follows : 


1.  Have  five  A/ C fly  directly  to  the  destination,  each  having  a load 


490 


5 98  k lbs.  Total  fuel  cost  = 5 CQ  (98  k). 


* W 

equal  to  w.  = — 

1 M . 
min 

2.  Have  three  A/C  refuel,  each  having  a load  equal  to  w*  = — = ^22.  = 

i Q 3 

163.3  k lbs.  The  cost  will  be  3 Cx  (163.3  k lbs.).  To  find  this, 
we  have  to  solve  problem  (P2)  for  wQ  = 163.3  k lbs. 

3.  Have  three  A/ C fly  directly  with  each  carrying  w*  = Z = 120  k lbs 

i o 

(i  = 1 »2 ,3) , and  a fourth  A/ C refuel  in  mid-air,  carrying  w*  = 

(W-3Z  )/l  = W - 3Z  = 130  k lbs.  The  total  cost  is  = 3 C (Z  ) 
u ° o o 
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+ ci  (w  “ 3Zq)  - 3 CQ  (120)  + Cj  (130).  The  quantity  (130)  can 
be  found  by  solving  problem  (P2)  for  wQ  = 130  k lbs. 

4.  Have  two  A/ C fly  directly  with  each  carrying  w*  = ZQ  = 120  k lbs 

( i — 1»2)»  and  another  two  refuel  in  mid- air,  each  carrying  w^  — (W— 
2Zq)/2  = 125  k lbs.  (i  = 3,4).  The  total  cost  is  = 2 CQ  (ZQ)  + 2C^ 
[(W-2Zq)/2]  = 2 Cq  (120)  + 2 (125).  The  quantity  (125)  can  be 

found  by  solving  problem  (P2)  for  wq  = 125.  Finally,  compare  the 
four  different  total  costs  and  select  the  minimum  one. 

Summary  of  the  Solution  Procedure  of  Problem  (P4) 

1.  Find  Mmin  according  to  eq.  (6-1). 

2.  Find  Q according  to  eq.  (6—2)  or  use  the  suggested  procedure  (found 
on  page  72)  if  Z^  is  not  known. 

3.  Set  = 0,  £r  = Q,  where 

&D  = the  number  of  A/C  flying  directly  to  the 
destination,  and 

% = the  number  of  A/°  that  need  to  refuel  in  mid- 
air 

4.  Let  *D  ■ *D  + 1.  *R  - *R  - 1 

5.  If  ^ = 0 go  to  step  8. 

6.  Check  the  feasibility  of  the  current  (j^D  £RR)  solution  as  follows: 

W “ ^DZo 

If  : < Z,  then  this  solution  is  feasible.  Record  this 

*R 

solution  and  go  to  step  4.  Otherwise  go  to  step  7. 

7.  Increase  £D  by  1,  i.e.  let  = JlD  + 1.  If  go  to  step 

8.  Otherwise  go  to  step  6. 


82 


8.  Compute  the  cost  of  each  feasible  solution  as  follows: 


(a)  For  the  nD)  solution  we  have: 


* W 
w.  = 


i M 


i = 1,2,  . . . M 


min 


min 


Total  cost  = M . * C — ) 

min  o 'M  ^ 
min 


(b)  For  the  (QR)  solution  we  have 


* 

w. 

l 


i = 1,2,  , . . Q 


Total  cost  = Q * 


Solve  problem  (P2)  with  wQ 


to  find  C^ 


(c) 


For  any  solution  of  the  form  (ID  £ R)  we  have: 

D R 


* 


w 


i 


Z 

o 


for  i = 1,2, 


» 


* 


w 


i 


W - 1Z 
D o 


for  i = £ + 1 , 


+ \ * 


W - £ Z 

Total  cost  = Ad  * CQ  (ZQ)  + *R  * Cj  ( . 

R 

W - LZ  w - V 

Solve  problem  (P2)  with  wQ=  — j to  find  Cj^  ( j 


9.  Compare  the  different  costs  and  pick  the  minimum  one. 


Figure  6-1  gives  a graphical  summary  of  the  solution  procedure 
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Record  this 
solution 


Yes 


Cd= 

^R  = ^R  _1 


feasible  solution 


Compute  the  cost  of  each  solution 
Compare  the  different  costs 
Pick  the  minimum 


Stop 


Figure  6-1:  Graphical  summary  of  the  solution  procedure  of 

problem  (P4) 


CHAPTER  VII 


PROBLEM  (P5) — MODEL  AND  SOLUTION 
1 A/C 

2 Refueling  Required 
(3  legs) 

Problem  Description 

This  chapter  deals  with  the  situation  where  the  transport  A/C  needs 
to  refuel  twice.  First  we  will  describe  the  problem  and  indicate  when 
this  need  arises.  Second,  we  will  provide  an  iterative  solution  pro- 
cedure for  problem  (P5),  which  depends  on  the  solution  of  problem  (P2). 

In  general,  there  are  two  cases  where  the  transport  A/C  needs  to 
refuel  twice.  The  first  case  involves  only  one  tanker  base  which  has 
two  tanker  A/C.  We  know  from  problem  (P2)  that  the  feasible  points  must 
lie  in  a region  whose  points  must  satisfy  eqs.  (4-1),  (4-2)  and  (4-3). 
However,  it  could  be  the  case  that  these  equations  do  not  simultaneously 
define  one  region.  Rather,  two  regions  are  formed:  one  region  where 

eqs.  (4-1)  and  (4-3)  are  satisfied,  and  another  where  eqs.  (4-2)  and  (4- 
3)  are  satisfied.  (See  Figures  (7—1)  and  (1—2).)  Therefore,  one  tanker 
can  fly  to  a point  in  the  first  region  to  meet  and  refuel  the  transport 
A/C.  It  goes  back  to  its  base,  while  the  transport  flies  into  the 
second  region.  There,  the  transport  meets  another  tanker  from  the  same 
base  to  refuel  for  the  second  time.  Afterwards,  it  heads  directly  for 
the  desired  destination.  Here,  aerial  refueling  is  done  twice,  using 
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’B’ 


Figure  7-1:  One  case  where  refueling  twice  is  necessary. 


36 


tankers  originating  from  the  same  base.  We  need  to  optimally  locate  the 
two  refueling  points  and  specify  the  initial  fuel  for  the  transport  and 
each  tanker. 

The  other  case  involves  two  tanker  bases  and  is  very  similar  to  the 
first  one.  (See  Figs.  (1—1)  and  (7—2).)  There,  the  tankers  originate 
from  two  different  bases.  In  our  treatment  here,  we  will  deal  with  the 
more  general  case,  namely  the  second  one  where  we  have  two  tanker  bases. 


Motivation  for  the  Solution  Procedure  of  (P5) 

We  will  solve  (P5)  in  an  iterative  way  where  each  iteration  con- 
sists of  two  steps  (2  subproblems).  At  each  step,  we  fix  one  refueling 
point  and  solve  for  the  other.  Then  we  use  the  results  of  the  first 
subproblem,  as  an  input  to  the  second,  to  resolve  for  the  original 
refueling  point.  This  process  is  repeated  and  if  convergence  can  be 
reached  then  the  solution  is  optimal.  Note  that  each  subproblem  is 
essentially  problem  (P2)  with  minor  changes. 


Let  "r 


1 


it  «t 

r2 


= the  first  refueling  point 
= the  second  refueling  point 
= distance  between  ”0”  and  "r 


1 


d; 

d2 


- the  distance  between  "r^"  and  "r ^ ' 
= the  distance  between  "r^”  and  ”B1" 
= d„ 


and 


d3 


= the  distance  between  "t2”  and  "D" 

= the  distance  between  ”r2”  adn  ”B2" 
(See  Fig.  (7-3)) 

Now  we  can  identify  the  following  two  subproblems: 
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Figure  7-2:  Another  case  where  refueling  twice  is  necessary 
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9o  — 9max^wo^ 


Figure  7-3:  Problem  (P5)  where  refueling  is  done  twice. 
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Subproblem  I 

Suppose  we  are  given  "r^".  Let  = the  maximum  amount  of  fuel  the 
transport  A/C  can  have  immediately  after  it  refuels  for  the  first  time. 

Then,  g^  g^x  (WQ)  + R^x  “ Fcf§niax  ^wo^»  wo>  dl^ 

- FC*(Hmax,  0,d3)  - FN*(0,d3)  (7-1) 

Letting  ”r^"  be  our  new  origin  base,  solve  for  , g' 

and  h'.  This  translates  to: 
o 

4 

(SP5I)  Minimize  £ f \ 

i=l  1 


(S.T.)  fj 

= FC  (g;,wo>dj) 

f 2 

= FN  (wo,d’) 

f3 

= FC*  (h^,0,d’) 

f4 

= FN*  (0,d’) 

0 

< gi  < *1 

0 

< ho  < Hmax 

4 

«;  + h;  > * fi 

d; 

< R [gT,w  ] 
I o 

di 

< R [G  (w  ) ,w  ] 

max  o o 

d; 

R*  (H  ,0) 

2 max 

(See  Fig.  (7-4)) 

Subproblem  II 

Suppose  we 

are  given  ’^2".  Let  g^  = the  amount  of  fuel  the  trans 

port  A/ C must  have  when  it  reaches  "r ^ order  for  it  to  be  able  to 
reach  D after  it  refuels  at  "r2"  (g^  could  be  equal  to  zero).  That  is 


gf  = Max  | FN(w  ,d ' ) + FC  (H  0,dl)  + FN  (0,dl)  - H 
+*  max  j j 


max 


, 0}  (7-2) 
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Figure  7-4:  Subproblem  I — solving  for  given  ”r^". 


max 
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Now  treating  "r2"  as  our  destination  and  requiring  that  the  transport 
A/ C ends  up  with  an  amount  of  fuel  equals  to  gf,  solve  for  "r ^ ” gQ  and 


hQ.  This  translates  to 


(SP5II) 


Minimize 


4 

Z f 
i-1  ] 


(S.T.)  fl  = FC  (g0,wQ,d1) 

f2  =FN(wo+gf,  d2 
f3  = FC*  (hQf  0,  d3) 
f4  = FN*  (0,  d3) 


0 < < g (w  ) 

o max  o 

0 < h < H 

o max 


§o  + ho  > s fi  + Sf 
i=l  r 

dl  < R [gmax(wo)’  Wo] 

d2  < R[Gmax(wo  + gf}»  V gf] 

do  < 4’  R*  [H  , 0] 

J 2 max’ 

(See  Fig.  (7-5)) 

Summary  of  the  Solution  Procedure  of  Problem  (P5) 

1.  Let  k = 0.  Choose  a very  small  number  e > 0 as  tolerance  for  the 
stopping  criterion. 

2.  Start  with  a point  that  satisfies  eqs.  (4-1)  and  (4-3).  Call  it  r°. 

3.  Let  k = 1. 

4.  Compute  g-j.  according  to  eqs.  (7-1). 

k— 1 k 

5.  Using  r , solve  subproblem  (SP5I)  to  find  r«,  g'  and  h* . 

i Zoo 

k 

6.  Using  x compute  g^  according  to  eq.  (7-2). 

k ]/ 

7.  Using  r0 , solve  subproblem  (SP5II)  to  find  r, , g , and  h . 

£ loo 
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Figure  7-5:  Subproblem  II — solving  for 


r^',  given  ”r2". 


93 


tc  Ir — 1 

8.  If  the  distance  between  and  r > e , let  k = k + 1 and  go  to 
step  (4).  Otherwise,  convergence  is  reached  and  the  solution 
is  r*.  r£,  gQ,  hQ,  g;  and  h\ 

What  can  we  say  about  the  solution  procedure  of  problem  (P5)?  Does 
it  converge?  If  it  does,  is  the  solution  a global  optimal  solution?  If 
we  can  not  show  that,  can  we  at  least  show  that  the  solution  is 
stable?  That  is,  if  the  process  happens  to  reach  the  global  optimal 
solution,  will  the  procedure  stop  there  or  will  it  drift  away  from  it? 
The  following  section  gives  us  some  facts  about  the  solution  procedure 
of  problem  (P5)  which  will  answer  the  questions  we  have  just  raised. 

Facts  About  Problem  (P5) 

1)  The  solution  procedure  of  each  subproblem,  namely  (SP5I)  and 
(SP5II),  converges  to  the  optimal  solution. 

Proof : 

See  the  proof  of  the  solution  procedure  of  problem  (P2)  because  each 
subproblem  is  essentially  problem  (P2). 

As  a consequence  of  this  fact,  and  in  order  to  help  us  talk 
about  the  convergence  of  the  solution  procedure  of  (P5),  we  restate 
steps  5 and  7 of  the  solution  procedure  as  follows: 

2)  a-  Given  point  r^  * , then  point  r^  is  the  best  refueling  point 

corresponding  to  it;  thus  it  is  the  optimal  solution  to  (SP5I) 
lc  Ic 

b-  Given  r^ , then  point  is  the  best  refueling  point  corresponding 

to  it;  thus  it  is  the  optimal  solution  to  (SP5II). 
k k-1 

3)  If  r^  ^ Tj,  , then  we  have  an  improvement. 
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Proof : 

We  need  to  show  that  the  total  fuel  consumption  at  iteration  k is  at 
most  equal  to  that  at  iteration  k— 1.  In  mathematical  terms,  we  need 
to  show  that 


+ + h'k  + g'k  < hk-1  + gk-1  + h'k_1  + s,k_1 


(7-3) 


Given  rk  1 and  the  fact  that  rk,  not  rk_1 , is  the  optimal  solution 
to  (SP5I),  we  have: 


h,k  + g,k  < h,k_1  + g,k-1 

o &o  o so 


(7-4) 


Now  add  the  fuel  consumption  at  point  rk  1 to  both  sides  of  (7-4)  to 
get: 


.•k  + + hk-‘  + g1;-1  < h'"-1  + + h*-‘ 


o + eo'1  <7-5> 


Now,  given  r£  and  the  fact  that  r*,  not  rk-1 , is  the  optimal 
solution  to  (SP5II),  we  have 


gk  + hk  < gk_1  + hk-1 
O O O O 


(7-6) 


Now,  add  the  fuel  consumption  at  point  rk  to  both  sides  of  (7-6)  to 
get 


gk  + hk  + g*k  + h*k  < gk-1+  hk_1  + g'k  + h'k 


(7-7) 


And  equation  (7-3)  follows  from  (7-5)  and  (7-7) 


Q.E ,0a 


4) 


a_  If  rk  1 = r1  , then  rk  = r 2 
, Tf  k * . k * 

b-  If  r2  = r2  , then  r^  = 
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Proof : 

This  follows  from  Facts  # 1 and  2.  Q.E.D. 

This  means,  if  the  procedure  reaches  any  of  the  global  optimal 
refueling  points,  it  will  immediately  give  us  the  other  one  and  the 
procedure  will  terminate  there  and  will  not  drift  away  from  it. 

Thus,  we  have  a stable  solution  procedure. 

Thus  far,  we  showed  that  the  procedure  will  at  least  converge  to 
the  local  optimal  solution,  and  if  it  reaches  the  global  optimal 
solution  it  will  stay  there.  But  there  is  evidence  it  might  converge  to 
the  global  optimal  solution.  A set  of  problems  (similar  to  those  of 
Table  4-2)  representing  a variety  of  possible  geographic  configuration 
and  complexities  were  solved.  All  converged  to  the  global  optimal 
solution.  Therefore,  we  conjecture  that  the  procedure  will  converge 
globally  for  real  world  problem. 


CHAPTER  VIII 


SUMMARY  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


Summary 


In  this  dissertation,  we  investigated  and  analyzed  the  general  A/C 
mid-air  refueling  problem  in  an  attempt  to  characterize  and  find  the 
optimal  solution.  It  was  found  that  eventhough  a mathematical  model  for 
the  general  problem  could  have  been  formulated,  it  would  have  provided 
little  insight  and  would  have  been  too  complicated  to  solve.  Instead, 
the  analysis  approach  used  in  this  dissertation  consisted  of  breaking 
the  general  problem  down  into  smaller  and  simpler  (sub) problems  that 
were  addressed  in  order  of  increasing  complexity.  Work  that  has  been 
done  regarding  this  problem  includes: 

1.  Deriving  the  fuel-related  functions,  which  are  summarized  in  Table 

2-1. 

2.  Solution  of  the  first  problem  (PI).  It  has  been  proven  that  the 
minimum  number  of  A/C  must  be  used  and  the  cargo  weight  should  be 
distributed  equally  among  those  A/C  used. 
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3.  Solution  of  the  second  problem  (P2).  S-convexity  of  the  problem  has 
been  shown  and  a solution  procedure  has  been  devised  that  utilizes 
the  relationship  between  the  various  decision  variables  at 
optimality.  Also,  a proof  of  the  convergence  of  this  procedure  is 
provided.  The  solution  of  (PI)  and  (P2)  was  essential  to  facilitate 
the  solution  of  the  other  problems. 

4.  Solutions  of  the  third  and  fourth  problems  (P3  and  P4).  They  are 

based  on  the  solutions  of  (PI)  and  (P2)  and  the  insight  obtained 

from  them.  The  main  contribution  here  was  to  find  the  optimal 
* 

weights  (w  fs)«  We  can  say  that  the  general  N A/C  problem  with  one 
refueling  has  been  solved  optimally. 

5.  Solution  of  the  fifth  problem  (P5).  We  proposed  an  iterative 
solution  procedure  where,  at  each  iteration,  one  refueling  point  is 
being  fixed  while  we  solve  for  the  other.  Then,  using  the  result  of 
the  first  subproblem  as  an  input  to  the  second,  we  resolve  for  the 
original  refueling  point.  Each  subproblem  is  essentially  problem 
(P2)  with  minor  changes.  We  showed  that  this  procedure  will  at 
least  converge  locally,  and  that  the  solution  is  stable.  Moreover, 
numerical  results  of  some  test  problems  support  the  conjecture  that 
this  procedure  will  converge  globally  for  real  world  problem.  The 
extension  of  this  procedure  to  the  multi-refueling  case  is  straight- 
forward. We  would  fix  all  refueling  points  but  one  which  we  solve 
for.  Then  we  repeat  this  process,  solving  for  a different  refueling 
point  at  each  time  until  convergence  is  reached.  A global 
convergence  proof  for  this  problem,  however,  is  beyond  the  scope  of 


this  dissertation 
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Problem  (P6)  and  (P7)  were  not  solved.  Note  that  if  each  w^  is 
decided  upon,  then  we  will  end  up  with  a combination  of  (PI),  (P2)  and 
(P5)  whose  solutions  are  known.  The  big  question  that  is  left  is:  what 

are  the  values  for  the  optimal  weights  (w.'s)?  As  a conjecture,  we 
believe  that  the  solution  to  problem  (P7)  might  be  similar  to  that  of 
(P4).  To  be  more  specific,  recall  the  definitions  of  ZQ  and  Z^  (Chap. 

V,  page  62),  the  definition  of  Mmin  and  Q (Chap.  VI,  page  77),  and  the 
definitions  of  and  &R  (Chap.  VI,  page  81).  Similarly,  define  Z2 , Q2 
and  l as  follows: 

Z2  = maximum  cargo  weight  one  transport  A/C  can  carry  from  ”0''  to 
"D"  by  refueling  twice  ( Z ^ > Z^). 

Q2  = minimum  number  of  transport  A/C  needed  to  transport  the  whole 
cargo  (W)  from  "0"  to  "D"  when  aerial  refueling  is  used  twice 
for  all  of  them  = ("w/Z^"]. 

£R2=  a certain  number  of  A/C  used  which  will  need  to  refuel  twice  in 
mid-air. 

We  conjecture  that  the  optimal  solution  will  be  one  of  the  following 
options: 

a)  The  D)  option  where  Mmi- n A/ C fly  directly,  each  with  w^  = 

«'V,' 

b)  The  (QR)  option  where  Q A/C  refuel  once  in  mid-air,  each 

* , 

with  w^  = W/Q. 

c)  The  (Q2  R2)  option  where  Q2  A/C  refuel  twice  in  mid-air,  each 

•k 

with  w^  = W/Q2« 

d)  The  (2.  D 2 R 2 R ) option(s)  where  2,  A/C  fly  directly,  2 _ A/C 

D K.  KZ  ZD  K. 

refuel  once  and  l ^ refuel  twice.  The  optimal  weight  for  the 
different  A/C  is: 
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w 


i 


Z 


for  i = 1,  2 . . . i 


o 


D 


* 


w 


i 


Z 


1 


for  i = £q+  1 


w 


* 


i 


for  1 3 £d+  &r  + 1 


Future  Research  Areas 


As  noted  above,  several  different  problems  are  still  open  for 
further  research.  Following  are  some  of  them: 

1.  Proving  the  global  convergence  of  the  solution  procedure  of  problem 


2.  Solving  problems  (P6)  and  (P7).  The  main  job  is  to  find  the 
relationship  between  the  different  weights  at  optimality; 

3.  Extending  the  research  to  the  case  where  there  are  more  .than  one 
type  of  transport  A/C  and/or  tanker  A/C; 

4.  Investigating  the  case  of  a forbidding  region  (i.e.,  region  where 
the  refueling  cannot  occur  and/or  the  A/ C cannot  fly,  such  as  a 
particular  country's  air  space). 


APPENDIX  A 


GLOSSARY  OF  NOTATIONS  AND  DEFINITIONS 


1) 

N 

= 

number  oE  available  A/ C 

2) 

n 

= 

transport  A/ C number  where  n = 1,2,...N 

3) 

EW 

= 

transport  A/ C Empty  Weight 

4) 

MTOW 

= 

Maximum  _Take-off  ^Weight  of  a transport  A/G 

5) 

F 

max 

= 

maximum  fuel  capacity  of  a transport  A/C 

u 

max 

= 

maximum  fuel  capacity  of  a tanker  A/G 

6) 

W 

= 

total  weight  of  the  cargo  to  be  transported 

7) 

wn 

= 

weight  of  the  portion  of  the  cargo  to  be 

loaded  in  the  nc  transport  A/C,  (a 

decision  variable)  where: 

N 

W = £ w 

, n 
n=l 

8) 

Sn 

= 

initial  fuel  of  the  n^  A/C  (a  decision 

variable) 

9) 

hk 

= 

f- 

initial  fuel  of  the  kL  tanker  where  k = 

1,2,...  ,K 

10) 

GWn 

= 

gross  weight  of  a transport  A/C,  where 

GWn  = EW  + gn  + wn 

GWmax 

= 

maximum  allowed  gross  weight  of  a transport 

A/C 

ID 

®max^wn^ 

= 

minimum  {Fm^  , MTOW  - w - Ew} 
max  ’ n 

Gmax^wn^ 

= 

minimum  {f  , GW  - w - Ew} 

max  max  n 1 

100 


101 


12)  MPF(GW) 

= distance  traveled  by  an  A/C  in  miles  per 
1 ,000  lbs  of  fuel  at  a given  gross  weight 
GW 

13)  R (g  ,w  ) 

n v°n*  a7 

ia  1- 

= Range  of  the  nLtl  transport  A/C  when  its 

14)  FC  (g  ,w  ,d) 

n &n’  n’ 

initial  fuel  is  gn  and  cargo  weight  is  wn 
= _Fuel  ^Consumed  by  the  ntn  transport  when  it 
flies  a distance  d,  its  initial  fuel  is  gn 
and  cargo  weight  is  w^ 

[Note:  d must  be  _<  Rn(gn>wn)  ] 

15)  FNn  K>d> 

= exact  amount  of  Fuel  Needed  to  fly  a 
distance  d when  the  cargo  weighs  wn 

16)  \ (W 

[d  must  be  < R (F  ,0)] 

= Range  of  the  kta  tanker  when  its  initial 

17)  FGk  (hk>vd) 

fuel  is  h,^  and  cargo  weight  is  w^  (usually 
wk-  = 0) 

= _Fuel  Consumed  by  the  k1-*1  tanker  when  it 
flies  a distance  d,  its  initial  fuel  is  hk 

18)  FN*  (wk>d) 

and  cargo  weight  is  w^ 

* 

Ld  must  be  < R^  (h^w^)] 

= exact  amount  of  Fuel  Needed  by  tanker  to 
fly  a distance  d when  its  cargo  weight  is 
wk  [d  must  be  £ r£  0^,0)] 

1 9 ) dod 

= shortest  distance  between  the  origin  base 
and  destination  measured  along  the  arc  of 
the  great  circle  connecting  the  two  points 

^0)  (0O,<*>O) 

= spherical  coordinates  of  the  origin  base 
(latitude  and  longitude) 
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( 0D » <t>o^ 

= spherical  coordinates  of  the  destination 

C 8B»<j>g) 

= spherical  coordinates  of  the  tanker  base 

( 0,  <j>) 

= any  point  on  the  sphere  (earth) 

1)  9,  4>) 

= spherical  distance  between  the  origin  base 
and  ( 9,  <j>) 

d 2 C 0»  <ji) 

= spherical  distance  between  the  destination 
and  ( 0,  <j>) 

d^(  0,  <j>) 

= spherical  distance  between  the  tanker  base 
and  ( 0,  <j))  . 

APPENDIX  B 


MEASURE  OF  DISTANCE 

Since  a suitably  small  portion  of  the  earth  can  be  considered  as  a 
plane,  the  shortest  distance  between  any  two  points  in  such  a region  is 
generally  approximated  by  Euclidean  distance.  However,  Litwhiler  and 
Aly  [22]  showed  that  the  Euclidean  assumption  for  even  moderately  large 
regions  (800  by  1000  miles)  can  lead  to  significant  error  in  objective 
value.  In  one  problem,  they  showed  that  using  the  Euclidean  assumption 
resulted  in  objective  error  of  about  18  percent  and  a location  error  of 
over  10,000  miles.  (See  Table  B-l.)  In  this  case  it  is  more 
appropriate  to  model  the  earth  as  a sphere  (even  though  it  is  not  a 
perfect  sphere,  but  rather  a spheroid). 

In  going  from  one  point  to  another  on  the  surface  of  a sphere 
(earth),  one  traverses  a curve,  not  a straight  line.  The  curve  having 
the  shortest  distance  between  two  points  on  the  surface  of  the  sphere  is 
the  great  circle.  The  great  circle  is  the  circle  in  E^  which  has  a 
center  at  the  center  of  the  sphere  and  is  drawn  on  the  surface  so  that 
its  radius  is  equal  to  the  radius  of  the  sphere.  (See  Figure  B-l.) 

Thus,  the  shortest  distance  between  two  points  on  a sphere  is  the  length 
of  the  shorter  great  circle  arc  connecting  them. 

Three  equivalent  methods  of  calculating  this  distance  have  been 
reported  in  the  literature.  All  the  methods  use  the  following  two 
facts:  (i)  the  great  circle  distance  between  two  points  on  a sphere  is 

directly  proportional  with  the  radius  of  the  sphere;  (ii)  the  distance 
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TABLE  B-l:  Error  in  distance  when  the  planar  assumption  is  used 


Angle 
between 
two  points 
(degrees) 

(1) 

Approximate 

distance 

using 

planar  assumption 
(nautical 
miles) 

(2) 

Exact 

distance 

(nautical 

miles) 

error  = 
(2)  - (1) 
(nautical 
miles) 

% 

error 
(2)  - (1) 

(2) 

15 

888.978 

891.524 

2.546 

0.286 

30 

1762.733 

1783.047 

20.315 

1.139 

45 

2606.303 

2674.571 

68.268 

2.552 

60 

3405.242 

3566.095 

160.853 

4.511 

75 

4145.871 

4457.619 

311.749 

6.994 

90 

4815.504 

5349.143 

533.638 

9.976 

105 

5402.678 

6240.667 

837.989 

13.428 

120 

5897.335 

7132.190 

1234.855 

17.314 

135 

6291.006 

8023.714 

1732.708 

21.595 

150 

6576.951 

8915.238 

2338.288 

26.228 

165 

6750.271 

9806.762 

3056.491 

31.167 

180 

6807.999 

10698.285 

3890.286 

36.364 
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Figure  B— 1 : The  Great  Circle. 
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between  two  points  on  a unit  sphere  is  identical  to  the  angle  (measured 
in  radians)  between  the  two  normals  of  the  sphere  at  these  points. 

Therefore,  distance  = R • a (B-l) 

where  R = the  radius  of  the  sphere  and 

a = the  angle  between  the  two  normals  at  these  two  points 
(measured  in  Radians).  (See  Figure  B-2.) 

The  three  methods  differ  only  in  the  way  they  find  the  angle  a,  but 
they  produce  the  same  result. 


Method  1 

Katz  and  Cooper  [20]  and  Litwhiler  and  Aly  [22]  considered  the 
distance  as 


d(X,Y)  = 2R  arcsin  (— UiL_lJLLL)  (b-2) 

Z R 

where 

X and  Y are  two  points  on  a sphere  with  radius  R and  centered  at 
the  origin, 

d(X,Y)  = the  distance  between  X and  Y measured  along  the  surface  of 
the  sphere,  and 
||*||  = the  Euclidean  norm. 


This  is  because 


Sin(|)= 


IX  - Y | 
2 


R 


_ I ]X  - Y| 
2R 


a . / X - Y U 

2 = arCSin  ( 2R~  } 


(B-3) 


(See  Figure  B-3.)  And  since  d(X,Y)  = Ra,  then  equation  (B-2)  follows 
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d(X,Y)  = Ro: 


Figure  B-2:  The  spherical  distance  between  two  points. 
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Figure  B-3:  Finding  the  angle  a using  method  1. 
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The  two  following  methods  use  a third  fact.  If  a is  the  angle 
between  any  two  nonzero  vectors  X and  Y,  then 


Now,  if  X and  Y are  on  a sphere,  then 
I I X | | = | | Y ||  =R 

Furthermore,  if  we  consider  a unit  sphere,  then 


R = 1 , thus  Cos  a = X • Y 

a = arccos  (X*Y)  (B-5) 

and  since  the  distance  = Ra , and  R = 1,  then  a itself  becomes  the 
measure  of  distance. 


Method  2 

By  representing  X and  Y in  spherical  coordinates  as  (R,  0 , <j>^)  and 

(R,92>^2^’  respectively,  instead  of  the  cartesian  coordinates 


and 

(y 

1 >72  >73 ) > 

then 

X1 

= 

R 

sin 

♦l 

cos  9 

x2 

= 

R 

sin 

*1 

sin  9 

x3 

= 

R 

cos 

*1 

7l 

= 

R 

sin 

*2 

cos  9 

72 

= 

R 

sin 

*2 

sin  9 

73 

= 

R 

cos 

* 2’ 

, and 

| |X| | = | | Y | | = R as  before,  so 
| |X| | • | |Y| | - R2,  and 
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X * Y = R [sin  sin  ^ cos  cos  02+  sin  ^l  sin  ^2  sin  91  sin  ®2  + 
cos  <t>^  cos  <f>2] 

2 

= R [sin  ^ sin  <j>2  (cos  9^  cos  9^  + sin  0^  sin  9^)  + cos  cos  ] 
By  recognizing  that 

cosCO^  - 9^)  = cos  6^  cos  0£  + sin  9^  sin  9^ , 
we  arrive  at 

cos  a = sin  <j>^  sin  ^ cos(92  - 0^)  + cos  ^ cos<}>2  ( B— 6 ) 

Notice  that  9 = the  geographical  longitude,  but  <(>  is  not  the  geogra- 
phical latitude.  In  fact,  cj>  = tt / 2 - latitude,  and  it  is  called  the  co- 
latitude. So,  if  one  would  rather  work  with  longitude  and  latitude 
instead,  then  equation  (B-6)  becomes 

cos  a = cos  <J>!  cos  (f>2  cos(92  - 9^)  + sin  <)>^  sin  <|>2  (8-7) 

where  = the  latitude  of  point  i.  Equation  (B— 7)  was  used  by  Drezner 
and  Wesolowsky  [12]  and  originally  reported  by  Donney  [9],  This  is  the 
formula  that  is  used  throughout  this  study. 

Method  3 

Lex  X = (x^ , X2 , x^)  be  a point  on  a unit  sphere.  We  can  represent 
X as  (1,0,<|>)  where: 

4>  = the  latitude  of  X 

9 = the  longitude  of  X,  then 

x1  = cos  <|>  cos  9 

x 2 = cos  <]>  cos  9 

X3 


= sin  <j> 
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Let  (a,b,c)  be  the  cartesian  coordinates  of  a point  Y on  the  surface  of 
a unit  sphere.  Then  the  distance  between  X and  Y,  using  equation  (B-5), 
is: 


d(X,Y)  = a = arccos  (X,Y) 

= arcos  (a  cos  cf>  cos  9 + b cos  <{,  sin  9 + c sin  <{>)  (B-8) 

Litwhiler  and  Aly  [22]  used  this  formula  to  characterize  the  shortest 
distance  between  any  two  points  on  the  sphere. 


APPENDIX  C 


SPHERICAL  CONVEXITY 

Sets  and  functions  defined  on  the  surface  of  a sphere  do  not 
satisfy  the  usual  convexity  definition  in  E . However,  this  definition 
has  been  extended  and  modified  [1,3,10,12]  for  the  case  of  spherical 
sets  and  functions.  Several  properties  and  results  of  spherical 
convexity  theory  have  been  developed  which  are,  to  a certain  extent, 
similar  to  the  convexity  theory  in  E^.  For  a complete  treatment  of  this 
subject,  see  Drezner  and  Wesolowsky  [12],  Aly  et  al.,  [1],  Drezner  [10] 
and  Ayken  [3].  Among  the  several  definitions,  results,  and  properties 
found  in  the  references  cited  above,  only  the  most  important  ones  are 
included  in  the  following  list. 

Definition  C-l  (Drezner  and  Wesolowsky  [12]) 

A set  of  points  on  a sphere  is  said  to  be  spherically  convex 
(s-convex  set)  if,  for  any  two  points  of  the  set,  the  whole  shorter 
great  circle  arc  connecting  them  is  included  in  the  set. 

Definition  C-2  (Drezner  and  Wesolowsky  [12]) 

A spherically  convex  combination  of  two  points  r^ , r^  is  defined  as 
the  point  r = p (r^r^.X)  that  lies  on  the  shorter  great  circle  arc 
between  r^  and  such  that  the  great  circle  distance  between  r^  and  r 
is  X * d(r1,r2)  for  ) 0 < X < 1. 

Definition  C-3  (Drezner  and  Wesolowsky  [12]) 

A function  f(r)  over  an  s-convex  set  is  said  to  be  a spherically 
convex  function  (s-convex  function)  if  for  every  0 < X < 1: 
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f(p(r^>r2>  S.  ^ + (1  - X)  f(r^)  for  any  element  and  of 

the  set. 

Theorem  C-l  (Drezner  and  Wesolowsky  [12]) 

Points  within  a circle  on  a unit  sphere  with  radius  < m/2  form  an 
s-convex  set. 

Theorem  C-2  (Drezner  and  Wesolowsky  [12]) 

The  distance  from  a given  point  r is  an  s-convex  function  within  a 
circle  of  radius  tt/2  and  center  r. 

Theorem  C-3  (Ayken  [3]) 

Every  minimum  of  an  s-convex  function  over  an  s-convex  set  is  a 
global  minimum. 

Theorem  C-4  (Ayken  [3]) 

If  an  s-convex  function  f(r)  over  an  s-convex  set  D has  two 
different  minima,  rp  r«  e D,  then  any  s-convex  combination  of  them  is 
also  a global  minimum. 

Theorem  C-5  (Ayken  [3]) 

An  s-strictly  convex  function  over  an  s-convex  set  D possesses  a 
unique  minima. 

Theorem  C-6  (Aly  et  al.  [1]) 

The  search  for  an  optimal  solution  to  the  spherical  Weber  problem, 
where  demand  points  are  not  located  entirely  on  a great  circle  arc,  can 
be  restricted  to  the  spherical  convex  hull  of  the  demand  points. 

Theorem  C-7  (Drezner  [10]) 

In  the  spherical  Weber  problem,  if  all  demand  points  are  located  on 
a great  circle,  so  is  the  optimal  solution. 
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Theorem  C-8  (Ayken  [3]) 

In  the  spherical  Weber  problem,  if  all  demand  points  are  included 
within  a spherical  circle  of  radius  n/4,  then  every  minima  is  a global 
minimum. 

Theorem  C-9  (Ayken  [3]) 

For  the  spherical  Weber  problem  with  m _>  3 and  3 > 0,  if  all 
demand  points  are  included  within  a spherical  disk  or  radius  tt/2  and  at 
least  three  demand  points  are  not  spherically  colinear,  then  the  problem 
possesses  a unique  minimum. 


APPENDIX  D 


SPECIAL  CONVEXITY  PROPERTIES  FOR  PROBLEM  (P2) 

Theorem  D-l 

The  function  FN  (w,d)  is  s-convex  over  a spherical  disc  of  radius 
<_  ir/4. 

Proof 

Let  r,  r^  and  r2  be  any  three  points  in  the  described  spherical 
disc.  Let  r3  = pCr^r^A)  for  some  A (0  < X < 1)  (See  definition  C-2) 
Also  let  d(r^)  = the  s-distance  between  r and  r^  for  i = 1,  2,  3.  Thus 
d(r^)  u/2  because  the  radius  of  the  spherical  disc  is  <_  tt/4 . Therefore 

d^)  is  a s-convex  function  -►  d[  p(r^  ,X)  ] < X d(r^)  + (1  - X)  dCr^) 
or  equivalently  d(r3)  <_  X dCr^  + (1  - X)  d(r2).  Since  FN  (w,d)  is  a 
convex  and  increasing  function  of  d,  then  applying  FN  to  both  sides  of 
the  above  inequality  results  in 

FN  [w,d(r3)]  < FN  [w,  XdCr^  + (1  - X)  d(r2)] 

<_  X FN  [w,d( r^ ) ] + (1  - X)  FN  [w,d(r?)] 

Thus  FN  [w,d]  is  s-convex  Q.E.D. 

The  following  theorem  is  an  extension  of  theorem  C-6  of  Aly  et  al. 

[1]  to  the  case  of  problem  (P2). 

Theorem  D-2 

The  search  for  the  optimal  refueling  point  of  problem  (P2)  can  be 
restricted  to  the  spherically  convex  hull  of  the  three  base  points  (the 
origin,  the  destination,  and  the  tanker  base). 
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Proof 


Let  V be  the  spherically  convex  hull  of  the  three  points  and  X be  a 
point  on  the  sphere  such  that  X £ V.  Let  d^  , d2  and  d^  be  the  spherical 
distance  between  X and  the  origin,  the  destination,  and  the  tanker  base, 
respectively. 

Aly  et  al.  [1]  proved  that  there  exists  a point  X'e  V that  domi- 
nates X in  the  spherical  distance  with  respect  to  the  three  base 
points.  That  is,  if  dj  , d^,  and  d^  are  the  distance  between  X*  and  the 
three  points,  then  dj^  < d^  for  i = 1,  2,  3. 

Since  fuel  consumption  increases  with  distance,  (i.e.,  FN  (..d^), 
FC(.,.,d^),  FN  (..d^)  and  FC  (.,.^3)  are  all  increasing  functions  of 
the  distance),  then 

the  objective  function  value  at  X' 

< the  objective  function  value  at  X. 


Q.E .D. 


The  following  lemma  will  be  used  later  to  arrive  at  a sufficient 
condition  which  guarantees  s-convexity  of  the  objective  function  of 
problem  (P2). 

Lemma  (D-l) 

In  order  to  prove  s-convexity  of  the  objective  function  of  problem 
(P2),  it  is  sufficient  to  show  that  at  any  feasible  point  (after 
transferring  the  North  Pole  to  it)  the  second  derivative  with  respect 
to  <)>  for  any  given  9 (direction)  is  non-negative. 


117 


Proof 

The  proof  will  be  given  in  three  (3)  parts:  Part  (a)  gives  a 

geometrical  interpretation  of  convexity  in  Euclidean  Space;  Part  (b) 
shows  how  to  apply  this  using  polar  coordinates;  Part  (c)  proves  the 
above  Lemma  using  analogy. 

(a)  The  definition  of  (regular)  convexity  in  Euclidean  Space  when 

cartesian  coordinates  are  used  states  that:  A.  function  f(x)  is  convex 

at  a point  x of  a convex  set  S if 

f [ Xx  + (1  - X)  x ] < X f(x)  + (1  - X)  f(x) 

for  each  Xs  (0,1)  and  for  each  xeS,  x * x (see  Reference  [5],  page 
109).  Note  that  Xx  + (1  - X)  x is  any  point  on  the  line  segment 
joining  x and  x.  (This  is  true  because  we  are  using  cartesian  coordi- 
nates, but  not  true  if  we  are  using  any  other  coordinate  system.)  This 
means  f has  to  be  convex  along  any  straight  line  passing  through  x,  no 
matter  what  its  direction  is. 

(b)  Thus,  proving  the  convexity  of  a function  at  a point  would  be 
easy  if  we  could  implement  this  geometric  insight. 

Now,  see  what  happens  when  we  transfer  the  origin  to  the  point 
under  consideration.  Using  polar  coordinates  (9,r),  a certain  line 
which  passes  through  the  (new)  origin  is  defined  by  9 = constant  = C, 
while  r is  allowed  to  vary.  Along  this  line  (or  direction),  the  func- 
tion is,  in  essence,  a function  of  one  variable,  namely  r.  Thus  the 
function  is  convex  at  the  origin  and  along  that  particular  line  if  we 


32f(c,r)  | 
3r2  r=0 


> 0 


have 
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Moreover  the  function  is  convex  at  the  origin  if 


92f(9,r)  | 

3r2  r=0 


> 0 


for  all  9 


(c)  Now  moving  back  to  the  case  where  the  function  is  spherical,  we 
can  prove  s-convexity  in  a similar  fashion. 

First,  note  that  the  great  circle  distance  is  the  shortest  distance 
between  any  two  points  on  the  surface  of  the  sphere,  and  it  is  analagous 
to  the  straight  line  distance  in  E . 

Second,  recall  the  following  two  definitions  (which  can  be  found, 
among  others,  in  Appendix  C): 


Def.  C-2  (Drezner  and  Wesolowsky  [12]) 

A spherically  convex  combination  of  two  points  r^ , r^  is  defined  as 
the  point  r=p  (r^,r2>X)  that  lies  on  the  shorter  great  circle  arc 
between  r^  and  r^  such  that  the  great  circle  distance  between 
r^  and  r is  X * d ( r ^ , r2 ) for  Xe  (0,1). 

Def.  C-3  (Drezner  and  Wesolowsky  [12]) 

A function  f(r)  over  an  s-convex  set  is  said  to  be  a spherically 
convex  function  (s-convex  function)  if  for  every  0 < X < 1: 

f [p(r1 ,r2,X)]  < X f (r^)  + (1  - X)  f(r2) 

for  every  element  r^  and  ^ of  the  set. 

This  implies  that  the  function  must  be  convex  along  any  great 
circle  track.  Thus,  in  order  to  prove  that  the  function  is  s-convex  at 
a point,  we  first  transfer  the  North  Pole  (NP)  to  that  point.  Note  that 
along  any  great  circle  arc  that  passes  through  the  (new)  NP,  we  have  0 = 
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constant  = C while  <}>  varies.  Thus  f becomes  a function  of  one  variable, 
namely  <t>  [see  Figure  ( D—  1 ) ] » 

So  f is  s-convex  at  NP  along  that  particular  great  circle  track  if: 

| >0 

9cJ>  <^=0 

And  in  order  for  the  function  to  be  s-convex  at  NP  we  must  have 

| > 0 for  all  0 Q.E.D. 

9<p  4?=0 

Condition  for  the  s-convexity  of  the  objective  function  of  problem  (P2) 

We  saw  from  Lemma  ( D— 1 ) that  the  objective  function  is  s-convex  if 
at  any  point  (after  transferring  the  North  Pole  to  it)  the  second  deri- 
vative with  respect  to  <j>  for  any  given  0 (direction)  is  non-negative. 
Therefore,  to  prove  the  s-convexity  of  the  objective  function,  we  need 
to  show  that  this  is  the  case  at  all  feasible  points  or  give  a condi- 
tions) under  which  it  is  true.  However,  a general  formula  for  the 
second  derivative  is  mathematically  untractable  and  does  not  allow  us  to 
draw  any  helpful  conclusion.  Therefore,  we  will  adopt  another  approach. 

We  want  to  construct  (manufacture)  a worst-case  situation  where  the 
second  derivative  at  some  point  would  attain  its  minimum  value  it  could 
ever  achieve.  If  it  is  non-negative  or  if  we  find  a condition  under 
which  it  is  non-negative,  then  we  are  done,  and  the  problem  is 
s-convex.  This  is  true  since,  by  construction,  all  other  cases  will 
produce  larger  values  of  the  second  derivative. 

This  can  be  doen  by  positioning  the  tanker  base  ("B"),  positioning 
the  refueling  point  and  by  selecting  the  direction  we  take  in  such  a way 
that  minimizes  the  second  derivative.  This  selection  is  based  on  some 
characteristics  inferred  when  looking  at  the  different  terms  of  which 


the  second  derivative  is  composed. 
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Figure  D-l : Transferring  the  North  Pole  to  a point  under 

consideration. 
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The  first  thing  we  found  was  that  the  value  of  the  second  deriva- 
tive is  reduced  considerably  if  the  refueling  point  is  positioned  over 
the  tanker  base  ("B").  Moreover,  if  "B”  is  within  a certain  distance 
(for  the  C-5A  A/C  this  distance  is  at  least  3,950  nautical  miles)  from 
the  destination  "D" , then  the  second  derivative  is  minimized  if  "B" 
itself  is  located  on  the  great  circle  arc  joing  the  origin  ("0")  and  the 
destination  "D”.  Note  that  the  position  of  the  minimum  point  simplifies 
the  expression  for  the  second  derivative  and  makes  it  easy  to  draw  some 
conclusions.  In  fact,  the  second  derivative  is  found  to  be  non- 
negative. Since  this  is  the  worst  possible  case,  given  that  "B"  is 
within  a certain  distance  from  "D,"  then  the  function  is  s-convex.  This 
is  exactly  the  condition  given  in  Theorem  (D-3).  However,  it  is  not 
easy  to  prove  s-convexity  if  the  feasible  points  are  greater  than  3,950 
nautical  miles  from  "D"  because  then  we  get  a much  more  complicated 
expression  for  the  second  derivative.  That  is  because  the  minimum  point 
no  longer  lies  on  the  great  circle  route.  In  that  case,  the  expression 
for  the  second  derivative  is  no  longer  simple  enough  to  draw  similar 
conclusions  from.  It  should  be  noted,  however,  that  3,900  nautical 
miles  is  roughly  the  maximum  range  for  a fully  loaded  C-5A. 

In  summary,  when  the  feasible  points  are  within  a certain  distance 
from  the  destination  ("D"),  the  minimum  value  of  the  second  derivative 
occurs  if  "B"  is  located  on  the  great  circle  joing  "0"  and  "D,”  and  if 
the  refueling  point  itself  is  positioned  on  top  of  "B."  That  is  why  the 
condition  for  convexity  seems  to  be  independent  of  the  distance  from  the 
refueling  points  to  either  "B"  or  ”0." 

Now,  we  state  this  condition  as  theorem  (D-3)  with  an  example  and  a 


proof 
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Theorem  (D-3) 

The  objective  function  of  problem  (P2)  is  s-convex  at  any  feasible 
point  that  satisfies  the  following  condition: 

4^  ’ i, + 2 v cd-u 

where 


(j>  = the  distance  (measured  in  radians)  between  that  point  and  the 

destination  ("D")  and 

a = aQ  + al  (EW  4-  wQ),  where 

EW  = A/C  Empty  Weight 

wQ  = cargo  weight 

R = Radius  of  earth,  and 

aQ  and  a^  are  the  parameters  of  the  MPF  function 
(see  equation  (2-1)  and  Table  2-1) 

That  is,  if  the  points  are  within  a specific  distance  ( 4>D ) from  the 
destination  then  the  function  is  s-convex  at  those  points,  where  ^ is 
the  solution  to  the  following  equation: 


2 

a 


-31R 


tan  4>d  + 2 4>d 


and 


D 


(D-2) 


Note  that  depends  on  the  cargo  weight  as  well  as  the  type  of  A./ C 

used. 

In  Table  (D-l),  we  present  values  of  <)>  which  correspond  to 
different  cargo  weights  (wQ)  for  the  C-5A  A/C. 
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TABLE  D-l:  The  distance  <|>  as  a function  of  cargo  weight  (wQ)  for  C5-A 


w0  (lbs) 

measured  in 

nautical  miles 

degrees 

radians 

0 

4,594 

77.3 

1.3497 

50  k 

4,481 

75.4 

1.3165 

100  k 

4,338 

73.0 

1.2746 

150  k 

4,160 

70.0 

1.2222 

200*  k 

3,950 

66.5 

1.1611 

*200  k lbs.  is  the  maximum  cargo  weight  for  the  C-5A  A/C  due  to 
structural  weight  limitations. 
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Note  that  the  maximum  value  of  we  would  hope  for  is  j/ 2 (or 


90°).  (See  theorems  (C-l)  and  (C-2)  of  Appendix  C.)  Before  we  derive 
condition  (D-l),  we  would  like  to  give  an  example  that  will  show  the 
implication  of  this  condition. 

Example 

Suppose  the  Origin  ("0")  and  Destination  ("D”)  are  90°  (or  Jl) 
apart.  Moreover,  the  cargo  weighs  200  k lbs.  We  would  like  to  know 
whethec  or  not  the  function  is  s-convex  over  all  the  feasible  region. 
Solution 

Since  the  two  bases  are  90°  apart  then 

Dqq  = 6,220  status  miles  = 5349.12  nautical  miles 

w = 200  k lbs 


o 


The  data  for  the  C-5A  are  as  follows: 


aQ  = 36.283 


ax  = -0.027 


MTOW  = 670  k lbs.,  GWm_  = 720  k lbs 

Ulcd  X. 


EW  = 320  k lbs 


F 


= 300  k lbs  of  fuel 


max 


= 150  k lbs  of  fuel 


R,  = R(w  . (w  ))  = 3,032  nautical  miles  (or  51°) 

1 o’  &max  o 


'max 


G (w_)  = Min  { F , GW  — EW  - w } 
max  v o'  l max’  max  oJ 


max 


max 


max 


= 200  k lbs  of  fuel 


R2  = R(w0,  Gmax  (wQ))  = 3,900  nautical  miles  (or  65.7°) 
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But  foe  the  C-5A  with  200  k lbs.  of  cargo,  we  have  ^ = 3950 
nautical  miles  = t56.5°,  and  since  R.2  (in  degrees)  = 65.7°  < <^,  then  the 
function  is  s-convex  over  all  the  feasible  region.  (See  Figure  D-2.) 
This  implies  that  for  any  A./C  we  can  have  a s-convex  function  even 
though  the  distance  between  the  origin  and  destination  is  as  large 
as  tt/2  (or  90°)  provided  that  <j^  < This  shows  that  what  is  impor- 

tant is  the  distance  from  the  feasible  points  to  the  destination,  not 
the  distance  from  the  origin  to  the  destination  (although  the  latter 
must  be  < tt/2  in  all  cases).  In  fact,  for  the  C-5A,  we  always  will  have 
a s-convex  function  over  the  feasible  region  regardless  of  the  cargo 
weight,  provided  that  it  is  not  more  than  200  k lbs.  and  the  distance 
between  the  origin  and  the  destination  is  not  more  than  tt/2.  (This  is 
derived  from  the  previous  table  and  the  corresponding  range  for  the 
different  weights.) 

Proof 

In  this  proof,  we  will  derive  condition  (D-l)  using  the  worst-case 
method.  We  will  construct  a situation  where  we  will  be  able  to  have  the 
minimum  value  ever  possible  of  the  second  derivative.  Then,  if  need  be, 
we  impose  any  condition  to  guarantee  its  non-negativity. 

How  is  this  done?  First,  we  study  the  expression  for  the  second 
derivative  and  infer  any  desired  characteristics  (or  requirements) 
needed  in:  (1)  the  position  of  "B" ; (2)  the  position  of  the  refueling 
points  and  (3)  the  direction  of  the  derivative  (or  9 ) in  order  to  make 
the  second  derivative  as  small  as  possible.  Second,  we  will  look  for 
the  points  and  the  direction  which  meet  these  requirements.  Finally,  we 
ascertain  whether  the  derivative  is  non-negative.  If  not,  we  find  out 


what  makes  it  non-negative. 
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Figure  D-2 : An  example  shows  the  implication  of  Theorem  D-3 
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4 

Lat  V(  9,  <j>)  = E fi  ( 9,  <}>) , then 
i=l 


av(  Q,  4>) 


derivative  of  V(  9,  <J>)  with  respect  to  <j> 
(evaluated  at  <j>  = 0)  ia  the  direction 
(represented  by  the  angle  9)  that  we  take 


i.e.,  the  derivative  at  the  north  pole  (where  ij>  - 0)  depends  on  the 


angle  9. 


3 V(M)  | 

2 ' 

94>  4>=0 

and  it  also  depends  on  9 


= the  second  derivative, 


3 'V(9,<|>)  | = _R_ 

9<J>2  4>=0  ST 

R 


„ a.  R 

2 2 1 
{cot  4>0  sin  (9  - 9q)  + cos~  (9  - 9Q)  g — } 


aiR 

H — ^2  {cot  sin2  (9  - 9^)  - cos'"  (9  - 9^)  g } 

/E2  2 

aiR 

+ — — {cot  <t>B  sin2  (9  - 9ft)  + cos'"  (9  - 9g)  “ g“} 


/E 


a,  R 


~ {cot  sin  (9  - 9b)  - cos2  (9  - 9g)  —■}  CD-3) 


where  a^  < 0 and 


= (a  + a]_  gQ)2  - 2 aL  R4> Q 

(0-4) 

= a2  + 2ax  R4>d 

(0-5) 

= (a1  + a.  h )2  - 2 a E# 
o 1 o 1 ts 

(D-6) 

= (a’)2  + 2a,Rfj> 

(D-7) 

Note  that,  in  this  case,  the  variables  are  the  g^s  and  4>  ^ f s , depending 
upon  the  location  of  the  point,  together  with  the  direction  9 we  take. 


(See  Fig.  (D-3).) 
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Figure  D-3:  The  values  of  0 , 0^'s  and  <j>^'s  for  a point  after 

transferring  the  North  Pole  to  it. 
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Also,  note  that  all  the  terms  of  the  second  derivative  are  greater 


than  or  equal  to  zero  except  the  following  two: 

2 

SL  R. 

[cos2(0  - 9 ) — - -]  = C0S2(  0 - 0 ),  and 

/E~  u al  z3/i  u 


L 


1 


R [ 2 ✓ ^ An  ^ 

— L COS  C9  * 9b)  — J = ~ 


a,  R 

V 


a.  R 


3/2  cos  ( e - 0B) 


because  a^  < 0. 


It  is  now  easy  to  find  the  desired  point  and  direction  by  noticing 
that  the  following  four  conditions  are  desired  in  the  points  and  direc- 
tion we  are  looking  for:  [See  equation  ( D— 3 ) » ] 


1,  have  cos2  (9  - 9^)  + 1 

2,  have  cos2  (0  - 9 ) 1 

jLJ 

3,  have  <j)  0 to  make  E.4-  and  E,  + (See  (D-6)  and  (D-7).) 

D J H 


l if  cot  *D  > -g- 
-a^R 

0 if  cot  *D  < — 

We  will  show,  in  incremental  fashion,  that  there  exist  points  and  a 

direction  that  have  all  these  characteristics.  The  first  two 

characteristics  are  satisfied  by  a point  along  the  great  circle  track 

joining  "0"  and  "B"  and  by  having  9 located  along  that  track  too  (see 

Figure  D-4).  That  is  9 = 0 + 18CP  and  8 = 8 or  0 = 0 . 

Moreover,  the  third  characteristic  will  also  be  satisfied  as  well 

if  we  let  <}>  = 0 i.e.,  choose  the  north  pole  to  be  at  "B.”  (See  Figure 

B 


4.  have  cos2  (9  - 9 ) + 


D-5 . ) Thus  , we  have 
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« 


B" 


Figure  D-4:  A point  along  the  Great  Circle  track 

joining  "0"  and  "B"  have  the  first  two 
desired  characteristics. 
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Figure  D-5: 


Point  "B"  (North  Pole)  has  the  first  three 
desired  characteristics. 
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£v(  Q,  <&  i _ aiR 

,3/2 


3.  R. 

a#  * + t3ln2(9" 


2 f 1 _ _1 , 

aiR  „3/2  „3/2  * 


a.K  „ 2 a 

= -ryr  + — {cot  4ksin*(  9 - ^ ) - cos“  ( 9 - — 

/e7 


i 


aiR 

2 


+ a^R 


{ I —l 

(a'  + ah  )3  a'3 

o 1 o o 


"alR  a2 

Note  that  the  condition  cot  > — — is  equivalent  to  — ■ — > 

(tan  <|>  -+-2  <+>  ) . This  represents  a region  around  ”D."  See  Figure  0-5. 

^ ° _alR  2 

For  the  region  where  cot  tj>^  > — — , we  will  have  cos  ( 9 - 9^)  - 1 

if  ”B"  (or  the  NF)  happens  to  be  on  the  line  joining  "0"  and  "D"  (see 

Figure  D-7).  In  other  words,  9^  = 9^  + 130°  and  0 = 0^  or  0 = 0^,  so 


32V(M)i  = flR  a’R 
2 ? I 

3 V <t>=0  Ex 


i R a.  R ~ i i 

JTl  - 372  + aiR  t 1 3 ■ T 1 

' E t'  (a’  + a.  h ) a* 

1 2 o 1 o o 


For  this  case,  the  minimum  is  achieved  by  setting  the  following  quan- 
tities as  follows: 


dl  = R^wo,go^  where 


<W  (wo>  lf  R(wo-fW)  5 V2 


30 


-a  V*2  ♦ 


-,  otherwise  (i.e.,  gQ  is  just  enough 
to  take  us  half  way) . 


Of  course,  d2  = DqD  - d^. 
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Figure  D-6:  A region  around  "D"  where  cot  <j>^  > 
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Figure  D-7 : Point  "B"  (North  Pole)  will  have  the  four  desired 

characteristics  if  it  is  located  along  the  line 
joining  "0"  and  "D.” 
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This  means  hQ  = E2;  that  Is,  supply  the  transport  A./ C with  all  the 
fuel  it  needs  for  the  second  leg. 


infill  , - aiE2  (. 


3(j)  <j>—0 


1 L_i 

( 2 . . , ,3/2  3 j 

(a  + 2a  d)  a 


+ axR 


1 T 
1 1 


3 ,3 


( a'  + a.  h ) a 


} 


o 1 o o 

Now  it  remains  to  show  that  the  above  expression  is  always  > 0, 


Note  that  = ^2* 


V* 


_^a+  -a  + 2ald2 

o a at 


. Writing  it  differently,  yields 


(a  + a^hQ)  = a + 2a^d2  • Therefore, 


, - a,R2  | 5 r 

d<t>  i>=0  (a  + a.h  ) 

10 


3} 


+ V2  1 1 3 ' -3} 

1 (a'  + a. h )J  a,J 


o 1 o' 
,3 


2 ,a3  - (a  + alho)3,  . _ ,a;3-  (ao  + alho)3, 

j **  ’a'3  (a-  '3l 


’ -V  l 3 3 

1 aJ(a  + a.h  T 
1 o 


+ a.h  )' 
0 o 1 o 


-2„,q  v 3a(a  + a.h  ) + a.2.  h^ 

3 v(.Ml|  = -a  R (-a  h ) [ — L_2l 

,.2  ^ I alR  ^ alV  L 3.  , .3  J 

3 cb  d>  =0  a (a  + a.h  ) 

1 o 


3a'  (a'  + a.  h ) + a'T  h^ 

+ a R (-a.h  ) [-2-2, t° i-2j 

1 lo  a’ 3(a'  + a.h  )3 

0 0 1 o 


(D-8) 


Now,  since  an  < 0 and  a = a'  + a.w  , then  a < a'  and  a + a^h 

* I o 1 o o J-  0 


< a'  + a,h  > 
o 2.  0 
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2,  , . .2 

a (a  + a,h  ) 

1  o 


a' 3(a'  + a.  h )3 
o o 1 o 


which  is  equivalent  to 


3a(a  + aihQ)  ^ 3ao(a;  + a^) 

~3  3 > 3 3 

a (a  + a,h  ) a'  (a'  + ah  ) 

1 o o o 1 o 


(0-9) 


and  since 


2 ,2 
a,  h 
1 o 


2  w2 
a,  h 
1 o 


3 3 

a (a  + a. h ) 


1 o 


3 3 

a (a'  + h ) 
o o 


(0-10) 


then  adding  (0-9)  and  (0-10)  yields 


j-3a(a  + aLhq)  + a*  h^  > 3*'Qtz'Q  + a^)  + ajh3 


a3(a  + a^hQ)3 


t’3 

o o 


a'  (a*  + a, h )' 
1 o 


(0-11) 


from  which  it  follows  after  looking  at  equation  (0-8)  that, 


92v(-!-^i  >o 

3<P  <j>=0 


Q.E.D 
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